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Abstract 


The reactions (1) H + O 2 = OH + O and (2) O + H 2 = OH + H are the most 
important elementary reactions in gas phase combustion. They are the main 
chain-branching reaction in the oxidation of H 2 and hydrocarbon fuels. 

In this study, rate coefficients of the reactions (1) and (2) have been 
measured over a wide range of composition, pressure, density and temperature 
behind the reflected shock waves. The experiments were performed using the 
shock tube- laser absorption spectroscopic technique to monitor OH radicals 
formed in the shock-heated H2/02/Ar mixtures. The OH radicals were detected 
using the Pi(5) line of (0,0) band of the A 2 Z + <— X 2 n transition of OH at 
310.023 nm (air). The data were analyzed with the aid of computer modeling. 
In the experiments great care was exercised to obtain high time resolution, 
linearity and signal- to-noise. 

The results are well represented by the Arrhenius expressions. The rate 
coefficient expression for reaction (1) obtained in this study is k\ = (7.13 ± 0.31) 
x 10 13 exp(-6957 ± 30 K/T) cm 3 mol- 1 s- 1 (1050 K < T < 2500 K) and a 
consensus expression for fcifrom a critical review of the most recent 
evaluations of (including our own) is k\ = 7.82 x 10 13 exp(-7105 K/T) 
cm 3 mol -1 s _1 (960 K < T < 5300 K). The rate coefficient expression for k 2 is 
given by k 2 = (1.88 ± 0.07) x 10 14 exp(-6897 ± 53 K/T) cn^mol-V 1 (1424 K < T 
< 2427 K). For k\, the temperature dependent A-factor and the correlation 
between the values of and the inverse reactant densities were not found. In 



the temperature range of this study, non-Arrhenius expression of k 2 which 
shows the upward curvature was not supported. 
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Chapter I. H + O2 = OH + O Reaction 


Introduction 

The reaction between hydrogen atoms and oxygen molecules 

H + 0 2 = OH + O A r H °298 = 70.2 kJ/mol ( 1 ) 

is the most important elementary reaction in gas phase combustion systems. 
It is the main chain branching reaction in the oxidation of H 2 and hydrocarbon 
fuels. About 80 % of the O 2 is removed from the main reaction zone of typical 
stoichiometric hydrocarbon-air flames at atmospheric pressure by this step 
[1], Because of the large reaction endothermidty, this reaction is likely to have 
a high activation energy, hence, this reaction is rate controlling even at high 
temperatures. Therefore, the induction delay and the flame speed are sensitive 
to the rate coefficient of this reaction. At high temperature and low pressure, 
reaction (1) produces O-atom and OH radical as the further chain branching 
and propagating carriers, respectively, according to the following reactions, 


0 + H 2 = OH + H 

A r H °298 = 7.8 kJ/mol 

(2) 

OH + H 2 = H 2 0 + H 

A r H ° 2 98 = ~ 62.8 kJ/mol 

(3) 


1 
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thus maintaining the chain branching process. Reaction (1) also determines 
the first and the second explosion limits of the H 2 /O 2 system which are 
primarily dependent upon the chain branching reaction. 

As reaction (1) plays a central role in the oxidation of H 2 and 
hydrocarbon fuels at high temperature, it has been the subject of numerous 
investigations and reviews. Baulch et al. [2] reviewed the rate coefficient data 
up to 1972 and recommended k\ = 2.2 x 10 14 exp(- 8450 K/T) cm 3 mol -1 s -1 over 
the temperature range 700 - 2500 K. Based on the molecularity and the 
activation energy of the reverse of reaction (1), they claimed that the 
activation energy would not be less than the endothermicity. In 1973 Schott 
[3] measured the exponential growth of chemiluminescence from CO + O = 
CO 2 * + h\ during the branched-chain ignition of H 2 /C 0 / 02 /Ar mixtures behind 
reflected shock waves and reported a rate coefficient expression with a 
temperature dependent A-factor, k\ = 1.22 x 10 17 T “ 0 907 exp(- 8369 K/T) 
cm 3 mol _1 s _1 in the temperature interval of 1250 to 2500 K. The reviews by 
Wamatz [1], Cohen and Westberg [4], and Cohen [5] essentially adopted 
Schott's results as their recommended rate coefficient expression. Thus 
Schott's rate coefficient expression was recognized as the most reliable one and 
was widely used in the modeling of hydrogen and hydrocarbon combustion 
systems. However, Schott’s rate coefficient data were questioned by Frank 
and Just [6] . They shock-heated H2/02/N20/Ar mixtures and initiated the 
reactions with O-atoms produced by thermal decomposition of N 2 O. The O- 
and H-atoms formed were monitored using Atomic Resonance Absorption 
Spectroscopy (ARAS). Reaction rate coefficient expressions of k\ and k$ were 
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obtained by simulating 0- and H-atom concentration profiles with a reaction 
mechanism of 10 elementary reactions. From their experimental results it 
was concluded that Schott's expression of k\ gave too low a value by a factor of 
more than 1.8 at temperature of 2000 K. Their results, however, agreed well 
with the recommended values of Baulch et al. [2] . A negative temperature 
dependence of the A-factor ofki was not found. Prompted by this experiment, 
the title reaction was reinvestigated by many authors using laser absorption 
spectroscopic or ARAS methods or both under various experimental 
conditions. Pirraglia et al. [7] reported k\ values obtained under pseudo-first- 
order conditions by the flash photolysis-shock tube technique that used the 
ARAS method to monitor H-atom concentrations in the temperature range 
from 962 to 1705 K. The Arrhenius rate coefficient expression reported was &i 
= 1.68 x 10 14 exp(- 8119 K/T) cm 3 mol" 1 s" 1 which was essentially the low 
temperature extension of Frank and Just's results. Masten et al. [8] measured 
ki values in the rich H^O^Ar system using laser absorption spectroscopy 
behind incident and reflected shock waves. OH radicals formed were monitored 
by utilizing the Ri(5) line in OH A 2 X + «- X 2 n transitions. The rate coefficient 
expression of k\ = 9.33 x 10 13 exp(— 7448 K/T) cm 3 mol 1 s 1 over the 
temperature range 1450 - 3370 K was reported. Yuan et al. [9] performed 
experiments using the same technique of Masten et al. but substituting the 
Pi(5) line for the Ri(5) line as a monitoring wavelength of OH radicals and 
presented a k\ expression with a temperature dependent A-factor similar to 
that of Schott, k 1 = 1.59 x 10 17 T - ° 927 exp(-8493 K/T) cu^mol^s" 1 in the 
temperature interval of 1050 - 2700 K. However the mixture compositions 
chosen were mostly stoichiometric or lean, from which the OH profiles 
generated were not sufficiently sensitive to derive k\ values with confidence. 
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Shin and Michael [10] studied reaction (1) using the laser photolysis-shock 
tube technique over the temperature range of 1103 - 2055 K. H-atoms were 
produced by photolyzing NH 3 or H 2 O with an exdmer laser and H-atom decay 
due to the reaction with a large amount of O 2 molecules present (pseudo-first- 
order condition) was monitored by ARAS. The rate coefficient expression 
reported was k\ = 6.93 x 10 13 exp(- 6917 K/T) cm 3 mol - 1 s -1 . Du and Hessler 
[11] also investigated reaction (1) together with the H 2 decomposition reaction 
at very high temperatures (2050 K < T < 5300 K) by using a laser flash 
absorption - shock tube technique to monitor the OH radicals formed behind 
incident shock waves. They combined their results with the analytical 
expressions of Shin and Michael, Masten et al. and Pirraglia et al. and 
presented a rate coefficient expression of &i = 9.76 x 10 13 exp(- 7474 K/T) 
cm 3 mol _ 1 s _1 over the wide temperature range of 960 - 5300 K. Recently Yu et 
al. [12] reexamined both data of Masten et al. and Yuan et al. by a rigorous 
optimization using two different model responses, the characteristic times(t25, 
t 5 o> ^ 75 ) and the time difference(t 75 - t 2 s)> and reported the source of the 
previous disagreement between their results to be the vibrational relaxation of 
02 - They extracted a new rate coefficient expression, k\ = 8.30 x 10 13 exp(- 
7253 K/T) cm 3 mol - 1 s _1 over the temperature range of 1336 - 3370 K. Yang et 
al. [13] also reinvestigated their previous work [14, 15] for k\ using a single 
beam laser absorption study of OH radical and then combined their data with 
that of Shin and Michael to obtain a new expression over the temperature 
range of 1100 - 3550 K, k x = 7.60 x 10 13 exp(- 7065 K/T) cnrhnoHs" 1 . 

Even though numerous experimental studies have been done by many 
investigators, there still remains significant discrepancies in k\ at high and low 
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temperatures, for example, factors of 2.5 and 1.5 discrepancies at 2500 K and 
1050 K, respectively. Recently Schott [16] compared the previously published 
k\ values with reactant densities at 2000 K and pointed out that there seemed 
to be a correlation between the derived k\ values and the inverse of reactant 
densities. The intermediate hydroperoxyl radical (HO 2 *) formed by reaction 
(1) can be stabilized to HO 2 by transferring excessive energy to collision 
partners at low temperature and high pressure (density). In this case, k\ 
would be density dependent. It is well known that HO 2 is thermodynamically 
uns table at temperatures above - 1300 K. If HO 2 radicals formed at high 
temperatures decompose immediately to the reactants (H + 02 ), not the 
products (OH + O), then it may be possible to obtain lower k\ values. 

Generally, reaction rate coefficients can be expressed in two different 
ways: two-parameter Arrhenius type, k = A exp(- 0/T); and three-parameter 
non- Arrhenius type, k = A T n exp(- 0/T). In the second expression, the A-factor 
is a function of temperature. The experimentally determined temperature 
dependence of the A-factor, n, was either 0 or in the range from - 0.927 to 2.0 
[3, 8, 9, 14]. Also, in theoretical studies, Miller [17, 18] calculated k\ and k~\ 
using quasiclassical trajectory (QCT) and quasiclassical trajectory - quantum 
mechanical threshold (QCT-QMT) methods based on the potential energy 
surface of Melius and Blint [19]. The calculation resulted in negative Arrhenius 
curvature and it was explained by the nonstatistical "recrossing effects at 
high temperatures. However, according to the recent theoretical calculation of 
Varandas et al. [20] for k\ using quasiclassical trajectory-internal energy 
qu an t um mechanical threshold (QCT-IEQMT) and an improved version of 
QCT-IEQMT methods which employed their fourth version of a double many 
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body expansion (DMBE IV) potential energy surface, the temperature 
dependence of the A-factor, n, was found to be zero in the temperature range of 
1000 - 3000 K. Hence, the curvature of the Arrhenius expression of k\ is not 
firmly established. 

In this study, a series of experiments over a wide range of composition, 
pressure, density, and temperature were performed behind reflected shock 
waves to address these problems mentioned above, i.e., (1) the values of k\, (2) 
the inverse density dependence of k\\ and (3) the temperature dependence of 
the A-factor. The reaction progress was followed by state-selective cw laser 
absorption spectroscopic measurements of OH radicals produced in shock 
heated H^O^Ar mixtures. The data were analyzed with the aid of computer 
modeling. In the experiment great care was taken to guarantee high signal-to- 
noise ratio, linearity, and time resolution. 


Experimental Section 

The experimental setup including optical, mechanical and electronic 
components is shown in Figure 1.1. 

A. Shock Tube System 

Experiments were performed in a rolled square stainless steel pressure 
driven shock tube with internal cross section of 63.5mm x 63.5mm. The length 


Li 



Figure 1. 1 Schematic diagram of the experimental setup. 
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of the driven section was 5.5 m and the driver section was 2.2 m. The driver 
section was separated from the driven section by an unscored dead-soft- 
temper al uminum diaphragm (Alufoil Products Co.) and pumped with a Welch 
Duo-Seal vacuum pump. A Varian zeolite-sieve was used to trap the oil 
backstreaming from the roughing pump. Pressure in the driver section was 
always maintained above 500 micron to further reduce effects of back 
streaming. The diaphragm was burst with a cross-shaped plunger. Mixtures 
of He and Ar served as the driver gas. 

The driven section was routinely pumped between experiments to 3 x 
10 -6 Torr by a Varian SD-700 roughing pump and a Varian Turbo-V60 pump 
equipped with a Varian SD-90 fore pump. The combined leak and outgassing 
rate was about 5 x 10“ 6 Torr/min. Driven section pressure was always 
maintained above 500 micron when pumped by the SD-90 fore pump in order 
to reduce backstreaming. The same type of sieve was used to trap the residual 
back flow of contaminants from the roughing pump to the shock tube. 

The gas handling system was constructed with high vacuum Varian 
valves (Model No. L6591-301) and seals. The vacuum pressure of the driven 
section and the gas handling system was measured by three Varian 0531 
thermocouple gauges coupled to Varian 801 vacuum readouts and a Vanan 
525 cold cathode gauge with a Varian 860A vacuum readout. Both the 
reactant pressures in the mixing tanks and the initial test gas pressures in the 
shock tube were measured by a Druck Model DPI-260 250 Torr pressure 
transducer (accuracy of 0.01 Torr) and DPI-260 5000 Torr pressure 
transducer (accuracy of 1 Torr, Druck, Inc.). 


H 
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Two 25.4 mm diameter Sl-UVA quartz windows were flush mounted on 
opposite sides of the shock tube's inner walls. Their centers were located 1.27 
cm from the reflective surface of the end plate. The shock velocities were 
measured by four 113A21 PCB Piezotronics pressure transducers mounted 
flush to the inner wall of the shock tube and powered by a PCB 482A05 
Piezotronics power supply. The distances between the pressure transducers 
were 45.72 cm, 30.48 cm, and 22.86 cm, respectively, downstream from the 
diaphragm. The last pressure transducer was positioned 1.27 cm from the end 
plate, center of the observation windows. 

Pressure transducer signals were fed into three two channel Nicolet 
digital oscilloscopes each set to 0.5 (is per point time resolution. Incident shock 
velocities measured between pressure transducer signals were fitted to a 
second order polynomial as a function of distance from the first pressure 
transducer. The shock velocity for the calculation of shock parameters at the 
reflected shock front was obtained by extrapolation of the fitted polynomial to 
the end plate. In this way the attenuation of shock velocity was automatically 
measured. Shock parameters were computed by standard methods [21] 
assuming full vibrational relaxation and no chemical reaction at the shock 
front and using NASA thermochemical data [22]. Later, the computed shock 
parameters were corrected for effects due to the interaction between the 
reflected shock front and the boundary layer built. 

B. Laser System 


The concentration of OH radicals during the course of reaction was 
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monitored using the Pi(5) line of (0,0) band of the A 2 £ + <- X 2 n transition at 
310.032 nm (air). A Coherent CR-699-21 ring-dye laser running with Kiton 
Red 620 dye (Exciton Chemical Co.) was pmnped by a Coherent Innova 200 
argon ion laser. The dye solution, Kiton Red 620 dye dissolved in methanol and 
ethylene glycol, was circulated by a Coherent 5920 Dye Circulator at 
pressures between 35 psi and 42 psi. The temperature of the dye solution was 
maintained at 5 °C by a Neslab CFT-25 Refrigerated Recirculator in order to 
achieve m aximum conversion efficiency. The concentration of the dye solution 
was controlled so that 90% of the argon ion laser beam power was absorbed in 
the dye jet. The argon ion laser is operated in single line (514.5 nm) and light 
regulation (7.5 ± 0.038 W) modes. The ring-dye laser was actively stabilized by 
adjusting the dye laser cavity length. An error signal generated by filtering out 
the amplitude information of a mode locked laser beam in a reference cavity 
was amplified and used to drive a high frequency piezoelectric mounted folding 
mirror and a scanning brewster plate to adjust the cavity length. Single 
frequency operation out of the CR-699-21 ring-dye laser was achieved by using 
a three-plate birefringent filter and an intracavity assembly (a thin etalon and 
a piezoelectric driven thick etalon). 

UV output (-310 nm) was generated via intracavity doubling of the 
fundamental beam using a Coherent Model 7500 frequency doubler (LiI 03 
crystal). The wavelength of the primary beam was monitored by a Burleigh 
WA-10 Wavemeter with a readout precision of ± 0.001 nm. The resulting 
single mode UV beam had a line width of 2 MHz and a typical output power of 
5 mW. The center wavelength of the OH absorption line was determined by 
passing part of the UV beam (see Detection System) through a burner 
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stabilized CH^Air flame. Laser wavelength was varied using the electronic 
controller until m axim um absorption within the flame was reached. At that 
time the 699-21 ring-dye laser was placed in a "Locked" mode. Frequency lock 
was then actively controlled by the laser electronics with no measurable 
deviation for the time between wavelength locking and shock arrival, i.e., 
approximately 30 seconds. The lasers, optical components, as well as the 
detection system described below were mounted on a pneumatically stabilized 
Newport MST series optical table. 

C. Detection System 

For signal detection a double beam scheme was employed. The UV 
beam was split into two beams using a UV beam splitter (Esco Products, Inc.). 
The first beam was directed through the shock tube at the center of the 
windows. If the expected absorbance for an experiment was small, the beam 
was reflected back through the shock tube at the same streamwise distance 
but vertically displaced from the incoming beam and onto the detector, i.e., 
double pass. For experiments with large expected absorbance, the beam was 
returned to the detector after being steered around the shock tube, i.e., single 
pass. An iris was set in front of the detector both to limit emission from the 
hot gases within the shock tube and to help establish the Time Zero 
measurement from the schlieren caused by shock passage. The second beam 
was split again and one part was directed to a reference detector and the other 
directed to a third detector through a CH^air flame. All three detectors were 
THORN EMI Model 9924QB Photomultiplier tubes (PMT) powered by THORN 
EMI Model PM28B power supplies. A customized five dynode configuration 
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was used for each PMT to guarantee optimum linearity for the light intensity 
of the UV laser beams. A high speed buffer/amplifier (National Semiconductor 
LH0033CG) isolated the anode from the coaxial signal cable. An overall 
electronic time constant of 0.2 ps was determined electronically for the entire 
PMT/electronic/cable system. Background light reduction was achieved by 
mounting a narrow band interference filter (Pomfret Research Optics, Inc., 
Model No. 20-3100-1, X 0 310 nm and AX at FWHM 0.1 nm) in front of each 
PMT. Also, during the course of experiments the room fights were extinguished. 
The reference signal, the probe signal, the difference (probe DC - reference 
AC), and the signal from the last pressure transducer were recorded by a four- 
channel Nicolet 4094C digital oscilloscope with Nicolet 4570 plug-in units and a 
Nicolet XF-44 disk recorder. The pressure signal served as the trigger source 
for the oscilloscope. Reaction progress was followed using the difference signal. 
The signal-to-noise ratio was typically 120-250. A laboratory PC (Dell 486P 
/33) running Asystant GPIB Program was used to control the Nicolet 4094C 
oscilloscope through an IOtech IEEE488 (Driver488, version 2.5) interface. 

Absolute pressure in the shock tube was measured using the voltage 
signal of the 113A21 PCB Piezotronics pressure transducer, coincidently 
mounted with the windows, and then converted to an absolute pressure using a 
calibration curve supplied with the transducer by the manufacturer (the 
calibration was validated against a dead weight gauge at the NASA Lewis 
Research Center Calibration Laboratory before use). The initial test gas 
pressure was added to this pressure. All shocks were initiated within one 
minute of filling the driver section with test gas. This actual pressure was used 
in the temperature correction for the boundary layer effect. 
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D. Test Gas Mixtures 

The test gas mixtures were prepared manometrically and allowed to 
stand for 48 hours before use. The maximum uncertainty of the final reactant 
concentrations was about 0.5% or less for each component. The punties of the 
gases stated were also confirmed by in-house gas chromatographic analysis 
using the flame ionization detector. They were: H 2 , 99.9995% (Linde Research 
Grade, THC as CH 4 < 0.3 ppm), 0 2 , 99.6% (Linde zero grade, THC as CH 4 < 
0.3 ppm), Ar, 99.996% (Linde zero grade, THC as CH 4 < 0.3 ppm). H 2 , 0 2 , and 
Ar were used without further purification. 


Results 

A. Experimental Conditions and Results 

Typical transmitted light intensity profiles at different experimental 
conditions are shown in Figures 1.2 (a) - (d). After an initiation - induction 
period, for min g small amounts of chain carriers H, 0, and OH, the light 
absorption rises rapidly due to essentially exponential growth of OH caused by 
the chain branching and chain propagating reactions. Depending upon 
experimental conditions the OH concentration reaches a maximu m and decays 
down later or reaches a partial equilibrium. 


In order to obtain the full information content of the absorption profiles 
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Figure 1.2(a) A typical experimental record of OH absorption at low 
temperature. Xh 2 = 0.02, Xo 2 = 0.005, Xa r = 0.975, T5 = 1234 K and P5 = 
0.945 atm. and P^max were chosen as simulation parameters. The 

schlieren signal corresponds to passage of the reflected shock front. The 
smooth line indicates the computed OH absorption profile using the Table 
1.3 reaction mechanism and the OH absorption coefficient obtained from 
self-calibration at maximum absorption (A ma3C ). The mismatch of the 
computed profiles after the peak do not affect determination of k\. See 
text for details. 
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Time(|Lis) 


Figure 1.2(b) A typical experimental record of OH absorption at low- 
middle temperature. Xh 2 = 0.02, Xo 2 = 0.005, X^r = 0.975, T5 = 1556 K and 
P 5 = 0.751 atm. NS^ and h max were chosen as simulation parameters. 
The schlieren signal corresponds to passage of the reflected shock front. 
The smooth line indicates the computed OH absorption profile using the 
Table 1.3 reaction mechanism and the OH absorption coefficient obtained 
from self-calibration at maximum absorption 
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Figure 1.2(c) A typical experimental record of OH absorption at high- 
middle temperature. Xh 2 = 0.01, Xo 2 = 0.001, Xat = 0.989, T5 = 1700 K and 
P5 = 3.203 atm. NS ma * and A max were chosen as simulation parameters. 
The schlieren signal corresponds to passage of the reflected shock front. 
The smooth line indicates the computed OH absorption profile using the 
Table 1.3 reaction mechanism and the OH absorption coefficient obtained 
from self-calibration at maximum absorption (A^o*). 
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0 25 50 75 100 125 

Time(fJs) 

Figure 1.2(d) Typical experimental and normalized slope records of OH 
absorption at high temperature. The upper panel shows an experimental 
profile of OH absorption at the condition of Xh 2 = 0.02, Xo 2 = 0.002, Xa r = 
0.978, T 5 = 2163 K and P 5 = 1.994 atm. The lower panel shows the 
normalized slope for the absorption profile in the upper panel. Smooth line 
represents the computed profiles using the Table 1.3 reaction mechanism 
and the OH absorption coefficient obtained from self-calibration at 
maximum absorption (A wcuc ). 
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while keeping the number of data points within reasonable bounds, five 
parameters were measured to characterize each experimental record [23], 
namely, A max = d - / / h)max', NS max = \{d(I / 1 0 ) / dt) max } I (s _1 ); three 
characteristic times, t 25 , tso, t 75 , at which the absorbed light intensity has 
reached 25%, 50%, 75% of A max , respectively. The value of NS max was 
obtained from each experimental profile in the following manner: the maxim u m 
slope ((dI/dt) mox ) was evaluated with the "Derivative" Standard Package 
Program for the Nicolet 4094 digital oscilloscope Gower panel of Figure 1.2 (d)) 
and normalized with I<>; then the absolute value was taken. 

The upper bound for the experimental temperature range was set by the 
appearance of OH absorption in a series of shocks in a 2.0% O 2 / 98.0% Ar 
mix ture. This mixture was chosen for its sensitivity to possible hydrocarbon 
contamination. At temperatures exceeding 2700K, a small but noticeable 
absorption was observed at long times. It is speculated that contaminants 
(possibly residual pump oil) swept off the shock tube walls or impurities in the 
test gas may be responsible. The lower temperature bound was set to 1050 K 
by the limit of the available observation time and the loss of signal due to the 
diminishing amounts of OH generated at low temperatures. Eight rich 
mixtures of equivalence ratio(<|>) of 2, 5, and 10 were used. Reflected shock 
temperatures ranged from 1050 to 2500 K and pressures ranged from 0.7 to 
4.0 atm. Mixture compositions were selected such that m a xim u m absorption 
was less than 0.55 for all mixtures. 

The temperature behind the reflected shock wave was corrected for the 
boundary layer effect with the measured absolute pressure (P 5 ) using the 
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adiabatic equation of state method described by Michael and Sutherland [24] . 
The corrected temperatures were always higher than the ideal shock 
temperatures by an average of 1.4%. A dependence of the temperature (T5) 
and the initial pressure (Pi) upon the corrected temperature was not observed. 

Experimental conditions and results of k\ are shown in Table 1.1. 

B. Computer Simulation 

Computer simulations were performed using the detailed reaction 
mechanism of Yuan et al. [9], with the following modifications. The rate 
coefficient expression for the OH + H2 reaction was taken from Ref. 25; and the 
H + OH + M reaction was added with the rate coefficient expression and third 
body efficiencies from Masten et al. [8]. The reaction mechanism and rate 
coefficient expressions are given in Table 1.2. Reverse reaction rate 
coefficients were automatically calculated from the principle of detailed 
balancing. NASA thermodynamic data [22] for the species were used in the 
calculation of shock parameters and equilibrium constants. The LSODE 
integrator [26] was utilized to solve a set of stiff differential equations 
describing chemical kinetics under assumed constant density conditions, 
calculated for the reflected shock waves [21]. 

The center line absorption coefficient of OH used to generate the OH 
absorption profile by computer simulation was computed according to the 
formula in the literature [9, 27, 28, 29, 30], namely: 
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Table 1.1 Experimental Conditions and Results of k\ a 


Ts 

Ps 

Amax 

NS;raax 

t25 

t50 

t75 

*i fc /10 12 



Xh 2 = 

=0.04, Xo 2 =0.01, XAr=0.95 (<t>=2) 


1052 

2.289 

0.334 

5246 

597 

618 

650 

0.092 

1074 

0.964 

0.303 

2210 

974 

1005 

1050 

0.104 

1086 

0.940 

0.294 

1953 

922 

985 

1053 

0.112 

1102 

0.960 

0.344 

2495 

802 

836 

861 

0.124 

1115 

2.248 

0.492 

9088 

371 

393 

411 

0.130 


Xh 2 =0.02, Xq 2 =0.005, Xat= 0.975 (<t>=2) 


1155 

0.957 

0.215 

999 

1188 

1274 

1343 

0.172 

1164 

0.892 

0.220 

1012 

1204 

1283 

1349 

0.179 

1170 

0.950 

0.232 

1134 

1111 

1183 

1245 

0.186 

1174 

0.945 

0.237 

1098 

1080 

1153 

1218 

0.181 

1228 

0.734 

0.227 

1057 

1079 

1153 

1217 

0.242 

1234 

0.945 

0.294 

1806 

768 

819 

863 

0.257 

1246 

0.877 

0.307 

1814 

753 

803 

848 

0.271 

1248 

0.756 

0.257 

1386 

886 

943 

989 

0.281 

1264 

0.773 

0.271 

1562 

856 

911 

961 

0.293 

1274 

0.781 

0.302 

1774 

761 

814 

860 

0.312 

1285 

0.794 

0.294 

1859 

738 

794 

843 

0.325 

1313 

0.936 

0.326 

2331 

572 

616 

653 

0.338 

1323 

0.800 

0.310 

2280 

649 

696 

737 

0.386 

1352 

0.766 

0.351 

2761 

581 

620 

658 

0.433 

1357 

0.771 

0.356 

2550 

572 

612 

651 

0.420 
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Table 1.1 Experimental Conditions and Results oik\ a (continued) 


t 5 

Ps 

Amax 

NS/720X 

t25 

t50 

t75 

k^/lO 12 



& 

X 

to 

II 

0.02, Xo 2 =0.005, Xat= 0.975 (<t>=2) 


1358 

0.772 

0.353 

2784 

590 

633 

670 

0.439 

1359 

0.774 

0.354 

2543 

571 

611 

643 

0.419 

1387 

0.778 

0.368 

2862 

500 

537 

576 

0.474 

1472 

0.859 

0.437 

4375 

355 

384 

412 

0.601 

1516 

0.699 

0.226 

2109 

426 

453 

498 

0.707 

1556 

0.751 

0.460 

5044 

302 

329 

354 

0.820 

1558 

0.733 

0.445 

4844 

316 

343 

369 

0.837 

1570 

0.737 

0.267 

3259 

287 

316 

347 

0.933 

1577 

0.962 

0.554 

7836 

240 

261 

280 

0.853 

1578 

0.788 

0.463 

5205 

283 

308 

333 

0.829 

1581 

0.761 

0.352 

4200 

257 

279 

303 

0.870 

1596 

0.852 

0.503 

6488 

248 

269 

290 

0.894 

1627 

0.738 

0.315 

4169 

227 

250 

275 

1.05 

1649 

0.733 

0.488 

6273 

248 

272 

296 

1.07 

1652 

0.769 

0.328 

4289 

214 

242 

262 

1.06 

1667 

0.744 

0.499 

6887 

232 

253 

276 

1.15 

1708 

0.771 

0.517 

7298 

208 

228 

248 

1.20 

1715 

0.797 

0.361 

5761 

198 

220 

241 

1.32 

1726 

0.784 

0.364 

5283 

196 

214 

235 

1.24 

1747 

0.795 

0.344 

5641 

186 

204 

223 

1.41 

1760 

1.034 

0.480 

9796 

132 

150 

163 

1.39 

1809 

0.846 

0.409 

7655 

146 

161 

177 

1.55 

1915 

0.876 

0.407 

9445 

118 

129 

141 

2.08 
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Table 1.1 Experimental Conditions and Results of k\ a (continued) 


Ts 

P 5 

A max 

NSmax 

t25 

t50 

t75 

fci fc /10 12 



Xh 2 =0.02, Xo 2 =0.005, Xat= 0.975 (<f>=2) 


1924 

0.735 

0.393 

6855 

130 

146 

164 

1.88 

1998 

1.018 

0.484 

12376 

90 

99 

111 

2.03 

2017 

0.757 

0.429 

9198 

114 

128 

143 

2.40 

2030 

0.854 

0.430 

10030 

101 

113 

125 

2.27 

2099 

0.772 

0.450 

9978 

97 

108 

121 

2.48 

2136 

0.851 

0.487 

12045 

81 

92 

101 

2.78 



Xh 2 =0.004, Xo 2 =0.001, Xat=0.995 (<J>=2) 


1525 

2.005 

0.210 

1206 

628 

703 

770 

0.734 

1639 

2.049 

0.270 

1821 

485 

550 

607 

1.00 

1692 

2.029 

0.284 

2115 

391 

432 

490 

1.17 

1934 

1.922 

0.341 

3430 

266 

293 

335 

2.07 

2044 

2.146 

0.438 

5582 

173 

199 

222 

2.52 

2068 

1.819 

0.390 

3841 

200 

236 

264 

2.32 

2082 

1.857 

0.375 

3836 

196 

224 

254 

2.35 

2132 

2.094 

0.439 

5835 

161 

185 

205 

2.83 

2256 

2.050 

0.439 

5956 

141 

163 

184 

3.10 

2409 

2.095 

0.453 

6978 

107 

127 

145 

3.67 



Xh 2 =0.002, Xo 2 =0.0005, Xat= 0.9975 (<t>=2) 


1527 

3.770 

0.156 

928 

573 

647 

704 

0.753 
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Table 1.1 Experimental Conditions and Results oik\ a (continued) 


T 5 Ps Amax NS max t 25 tso t75 ^1 6 /10 12 


Xh 2 =0.002, Xo 2 =0.0005, Xat=0.9975 (<)>=2) 


1593 3.704 

1658 3.654 

1716 3.792 

1725 3.905 

1763 3.843 

1888 3.880 
1977 4.050 

1999 3.697 

2001 4.088 

2092 3.681 

2211 3.715 


1243 0.798 

1246 0.795 

1251 0.792 

1252 0.800 

1254 0.813 

1301 0.754 

1315 0.766 

1330 0.784 

1333 0.793 


0.143 

961 

0.168 

1164 

0.217 

1580 

0.233 

1783 

0.195 

1595 

0.242 

2194 

0.271 

2862 

0.224 

2102 

0.276 

3270 

0.257 

2938 

0.270 

3433 


0.084 

576 

0.086 

630 

0.091 

662 

0.084 

595 

0.089 

676 

0.101 

818 

0.112 

961 

0.108 

1009 

0.108 

1025 


481 

534 

433 

476 

352 

397 

319 

360 

312 

356 

239 

270 

176 

203 

208 

237 

179 

203 

166 

193 

137 

158 


804 

863 

781 

822 

733 

784 

760 

807 

753 

799 

597 

643 

539 

581 

542 

578 

540 

578 


592 0.912 

515 1.09 

445 1.24 

406 1.27 

390 1.41 

302 1.77 

229 2.09 

264 2.08 

229 2.35 

217 2.66 

178 3.14 


907 0.259 

866 0.272 

820 0.276 

848 0.267 

838 0.282 

680 0.337 

614 0.359 

606 0.387 

610 0.391 


Xh 2 =0.05, Xq 2 =0.005, Xat=0.945 (<|>=5) 
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Table 1.1 Experimental Conditions and Results of k\ a (continued) 


T 5 Ps ^max NS/nax t25 t50 t75 A>1 6 /10 12 


Xh 2 =0.05, Xo 2 =0.005, Xat =0.945 (<>=5) 


1336 

0.797 

1336 

0.794 

1337 

0.796 

1516 

0.768 

1538 

0.794 

1543 

0.744 

1546 

0.799 

1549 

0.794 

1565 

0.810 

1620 

0.759 

1637 

0.760 

1675 

0.749 

1684 

0.758 

1706 

0.739 

1743 

0.749 

1754 

0.702 

1805 

0.756 

1897 

0.761 

1939 

0.728 

2014 

0.721 

2015 

0.746 

2038 

0.849 

2041 

0.732 


0.110 

1061 

0.107 

964 

0.114 

1098 

0.168 

2369 

0.152 

2252 

0.164 

2181 

0.152 

2333 

0.178 

2893 

0.198 

3219 

0.196 

3095 

0.171 

3094 

0.160 

2969 

0.199 

3847 

0.160 

2941 

0.231 

4392 

0.230 

4009 

0.264 

5378 

0.293 

7001 

0.260 

6414 

0.265 

6173 

0.290 

7531 

0.311 

9328 

0.298 

7081 


532 

574 

526 

559 

527 

561 

252 

274 

260 

280 

245 

269 

248 

269 

236 

257 

213 

235 

201 

220 

199 

216 

184 

200 

173 

191 

174 

190 

154 

168 

163 

177 

128 

140 

109 

121 

107 

118 

94 

105 

89 

102 

79 

86 

87 

98 


602 

0.368 

586 

0.392 

593 

0.395 

294 

0.722 

296 

0.755 

291 

0.725 

286 

0.777 

273 

0.822 

250 

0.833 

236 

0.910 

231 

1.03 

214 

1.13 

205 

1.17 

203 

1.18 

179 

1.26 

191 

1.26 

155 

1.43 

133 

1.80 

130 

2.02 

117 

2.11 

112 

2.25 

95 

2.34 

108 

2.18 
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Table 1.1 Experimental Conditions and Results ofk\ a (continued) 


Ts 

Ps 

Amax 

NS,™* 

t25 

t50 

t75 

fci fc /10 12 



Xh 2 =0.05, Xq 2 =0.005, Xat =0.945 (< t >=5) 


2125 

0.721 

0.324 

8794 

75 

88 

95 

2.69 

2151 

0.724 

0.337 

8907 

79 

84 

100 

2.69 

2188 

0.775 

0.343 

9794 

65 

74 

83 

2.78 

2376 

0.756 

0.383 

12817 

50 

58 

66 

3.69 

2414 

0.725 

0.366 

11792 

49 

57 

65 

3.84 



Xh 2 =0.02, Xo 2 =0.002, Xat =0.978 (<>=5) 


1509 

1.629 

0.098 

1259 

307 

338 

361 

0.736 

1522 

1.754 

0.107 

1430 

278 

313 

334 

0.760 

1531 

1.673 

0.108 

1513 

294 

321 

340 

0.802 

1578 

1.682 

0.126 

1641 

240 

264 

285 

0.808 

1667 

1.725 

0.154 

2477 

188 

209 

224 

1.06 

1746 

1.744 

0.170 

3527 

149 

165 

181 

1.44 

1819 

1.757 

0.195 

4428 

127 

144 

155 

1.69 

1850 

1.759 

0.184 

3775 

121 

136 

148 

1.60 

1884 

1.893 

0.223 

5429 

104 

116 

128 

1.81 

1948 

1.834 

0.215 

5174 

93 

106 

117 

1.98 

1989 

1.579 

0.201 

4373 

101 

113 

126 

2.16 

1997 

1.838 

0.229 

5587 

86 

97 

107 

2.11 

2026 

1.786 

0.244 

6289 

83 

95 

104 

2.34 

2029 

1.613 

0.219 

5124 

92 

106 

115 

2.36 

2049 

1.578 

0.197 

4731 

91 

103 

114 

2.50 
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Table 1.1 Experimental Conditions and Results of k\ a (continued) 


Ts 

Ps 

Amax 


t25 

t50 

t75 

ki b /10 12 



X H 2 =I 

0.02, Xo 2 =0.002, XAr=0.978 (<t>=5) 


2094 

1.816 

0.271 

7590 

70 

81 

90 

2.65 

2134 

1.783 

0.294 

8239 

67 

78 

87 

2.80 

2163 

1.994 

0.323 

10388 

56 

64 

72 

2.95 

2165 

1.981 

0.294 

9250 

59 

67 

75 

2.91 

2184 

1.780 

0.289 

8400 

60 

71 

80 

3.02 

2282 

1.849 

0.265 

8422 

51 

59 

67 

3.41 

2293 

1.852 

0.303 

9469 

50 

58 

65 

3.38 

2345 

1.915 

0.358 

12328 

49 

52 

60 

3.70 

2379 

2.024 

0.372 

13641 

40 

47 

54 

3.81 

2501 

2.013 

0.356 

14388 

34 

40 

46 

4.46 



ft 

to 

II 

0.01, Xo 2 =0.001, Xat=0.989 (0=5) 


1521 

3.263 

0.090 

1077 

298 

328 

352 

0.724 

1533 

3.303 

0.075 

976 

291 

323 

345 

0.760 

1557 

3.382 

0.088 

1281 

255 

280 

298 

0.850 

1579 

3.054 

0.104 

1359 

258 

286 

308 

0.888 

1583 

3.510 

0.106 

1455 

220 

246 

270 

0.818 

1609 

3.387 

0.091 

1203 

216 

241 

261 

0.833 

1611 

3.720 

0.109 

1672 

202 

224 

239 

0.881 

1647 

3.342 

0.100 

1597 

211 

231 

247 

1.05 

1664 

3.868 

0.123 

2241 

165 

185 

201 

1.07 

1700 

3.203 

0.117 

2004 

180 

201 

217 

1.25 


ii 
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Table 1.1 Experimental Conditions and Results of k\ a (continued) 


T 5 Ps ^max NS max tfcs tso t 75 ^l 6 /10 12 


Xh 2 =0.01, Xo 2 =0.001, Xat=0.989 (<>=5) 


1714 

3.237 

1792 

3.275 

1809 

3.331 

1841 

3.811 

1846 

3.364 

1895 

3.285 

1913 

3.378 

1945 

3.117 

1976 

3.280 

2057 

3.329 

2104 

3.439 

2107 

3.403 

2142 

3.444 

2155 

3.383 

2213 

3.370 

2305 

3.442 

2355 

3.478 

2379 

3.313 

2413 

3.384 

1514 

0.722 


0.120 

1957 

0.144 

2958 

0.148 

2891 

0.157 

3219 

0.150 

3032 

0.179 

3492 

0.183 

3716 

0.195 

3843 

0.189 

3928 

0.195 

5030 

0.253 

6653 

0.251 

6045 

0.208 

5561 

0.223 

5730 

0.224 

6294 

0.255 

8044 

0.229 

7196 

0.249 

7732 

0.262 

8553 


176 

197 

147 

165 

139 

153 

108 

122 

123 

138 

114 

130 

104 

118 

110 

126 

99 

111 

82 

93 

71 

82 

73 

84 

65 

77 

65 

77 

56 

67 

50 

59 

45 

53 

47 

56 

41 

49 


214 

1.19 

178 

1.59 

168 

1.53 

136 

1.46 

153 

1.61 

143 

1.72 

131 

1.76 

140 

1.91 

122 

1.97 

103 

2.55 

93 

2.65 

95 

2.49 

86 

2.78 

87 

2.75 

76 

3.13 

66 

3.65 

60 

3.73 

64 

3.91 

57 

4.10 


Xh 2 =0.10, Xo 2 =0.005, Xat=0.895 (0=10) 

0.075 1131 269 288 307 0.746 


Table 1.1 Experimental Conditions and Results of k\ a (continued) 


T 5 P 5 Amax NS mC x t25 ^50 t 75 ^i^/lO 1 ^ 


Xh 2 =0.10, Xo 2 =0.005, Xai -=0.895 (<fr=10) 


1569 

0.804 

0.096 

1714 

1575 

0.809 

0.093 

1783 

1579 

0.718 

0.091 

1538 

1616 

0.796 

0.086 

1755 

1698 

0.812 

0.115 

2865 

1715 

0.827 

0.117 

2960 

1782 

0.978 

0.136 

3988 

1801 

0.784 

0.132 

3341 

1869 

0.851 

0.169 

4725 

1883 

0.850 

0.148 

4527 

2042 

0.836 

0.143 

4876 

2081 

0.861 

0.157 

5804 

2407 

0.897 

0.206 

9744 

2481 

0.902 

0.211 

9930 

2494 

0.856 

0.215 

9787 


207 

229 

244 

0.845 

200 

217 

231 

0.904 

207 

226 

245 

0.900 

176 

193 

205 

1.01 

142 

155 

167 

1.30 

138 

152 

162 

1.32 

94 

104 

113 

1.39 

114 

127 

138 

1.52 

94 

104 

112 

1.67 

91 

102 

110 

1.82 

65 

74 

81 

2.41 

60 

68 

74 

2.63 

33 

39 

44 

4.16 

31 

36 

41 

4.31 

31 

37 

42 

4.42 


a Units are K for T 5 , atmosphere for P 5 , s 1 for NSn^*, M-S for t25, tso, and t 75 
and cm 3 mol - 1 s -1 for k\. 


b *l(NS nax) 
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Table 1.2 Reaction Mechanism for k\ a 



Reaction 

A 

n 

0 

Reference 

1. 

h + o 2 =oh + o 

7.13E+13 

0.0 

6957 

This study 

2. 

0 + H 2 = OH + H 

1.87E+14 

0.0 

6854 

31 

3. 

OH + H 2 = H 2 0 + H 

2.14E+08 

1.52 

1736 

25 

4. 

0 + H 2 0 = OH + OH 

4.51E+04 

2.70 

7323 

32 

5. 

o + o + m = o 2 + m 

1.00E+17 

-1.0 

0 

33 


Ar=1.0, H 2 =2.9, 0 2 =1.2, H 2 0= 

=18.5 




6. 

h + h + m = h 2 + m 

6.40E+17 

-1.0 

0 

33 


Ar=1.0, H 2 =4.0, H 2 O=12.0, H: 

=26.0 




7. 

H + 0 + M = OH + M 
Ar=1.0, H 2 O=5.0 

6.20E+16 

-0.6 

0 

8 

8. 

H + OH + M = H 2 0 + M 
Ar=1.0, H 2 =2.5, H 2 0=16.25 

8.40E+21 

-2.0 

0 

33 

9. 

h + o 2 +m = ho 2 + m 

7.00E+17 

-0.8 

0 

33 


Ar=1.0, H 2 =3.33, 0 2 =1.33, H 2 0=21.3 




10. 

ho 2 + H = OH + OH 

2.20E+14 

0.0 

710 

34 

11. 

ho 2 + h = h 2 + 0 2 

2.50E+13 

0.0 

350 

1 

12. 

ho 2 + H = H 2 0 + 0 

5.00E+12 

0.0 

710 

34 

13. 

H0 2 + 0 = 0 2 + OH 

2.00E+13 

0.0 

0 

1 

14. 

ho 2 + oh = h 2 o + 0 2 

2.00E+13 

0.0 

0 

1 

15. 

ho 2 + ho 2 = h 2 o 2 + 0 2 

1.06E+11 

0.0 

-855 

35 

16. 

h 2 o 2 + m = oh + oh + m 

Ar=0.67, O 2 =0.78, H 2 O=6.0 

1.20E+17 

0.0 

22900 

33 

17. 

h 2 o 2 + h = ho 2 + h 2 

1.70E+12 

0.0 

1900 

1 

18. 

H 2 0 2 + H = H 2 0 + OH 

1.00E+13 

0.0 

1805 

1 

19. 

h 2 o 2 + o = ho 2 + oh 

2.80E+13 

0.0 

3225 

35 

20. 

h 2 o 2 + oh = h 2 o + ho 2 

7.00E+12 

0.0 

720 

1 


a Rate coefficients are in the form k = A T 11 exp(- 0 / T). 
Units are cm 3 , K, mol, and s. 
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£cob = S x g(o) 0 ) 

= ^4^fiu(l - e-WksT)^), (4 ) 

m e c z 11 

where is the spectral line absorption coefficient at 0 ) o (s _1 ), S is the 
integrated absorption, g(co 0 ) is the value of the shape function (Voigt profile) at 
a>o, e and m e are the electronic charge and mass, respectively, c is the speed of 
light, n and n" are the total and ground state number densities of OH, 
respectively, kB and h are the Boltzmann and Planck constants, respectively, 
and fiu is the line absorption oscillator strength. The oscillator strength, fi u , 
was calculated from the following expression [36, 37, 38] : 

f, _ fuV'Sj'J-'Tj'j" (C) 

flu ” 2J" + 1 * (5) 

where f v ' v « is the v' <—v" band oscillator strength, Sj-j- is the J' 4- J" rotational 
line intensity factor (Honl-London factor), and Tj-j- is the correction factor for 
vibration-rotation interaction. For the Pi(5) line of the (0,0) band of OH A 2 Z + 
X 2 n transition employed in this study, f v ' v " = 1.1 x 10 -3 [39], Sj'j" = 6.13 [40], 
and Tj'j" = 0.978 [41] were used. The fraction of molecules in the ground state 
(n'Vn) was computed from the Boltzmann equation with the local thermal 
equilibrium (LTE) assumption, 

n " _ (2J" + 1) exp [- (EUEj) / k B T] 

n q 

where E„ and Ej are the vibrational and rotational energies, respectively, and 
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Q is the total internal partition function (Q„j x Q e i). The rovibrational 
partition function, Q^j was calculated using higher approximations for the 
rotational and vibrational energy levels [42], i.e., including terms for the 
anharmonicity of the potential energy curve, the influence of centrifugal force, 
and the fact that a molecule is a vibrating rotator. The reference energy level 
was defined as the u=0 state. In the rotational and vibrational energy 
calculation, the vibrational frequency (co e ), the anharmonicity correction (o)eX e ), 
the rotational-energy constant (B e ) at the equilibrium nuclear distance and the 
term for vibrating rotator (a) were taken from the literature [43]. For the 
electronic partition function, Q e i = 4 was assigned because of the Lambda- 
doubling of the ground n state of OH due to the coupling of spin-orbit and spin- 
intemuclear axis rotation [43] . 

In generating OH absorption profiles, because the information for the 
collision broadening for the Pi(5) fine is not known, the Voigt profile (g((0 o )) was 
not calculated. Instead, a self-calibration scheme was employed, i.e., the value 
of the shape function (g(cD 0 )), hence, e(OH), was systematically varied until 
satisfactory matches to the experimental observables were obtained. The 
absorption coefficients obtained in this way are shown in Figure 1.3. 

The local logarithmic response sensitivities defined below [44] were 
computed for the four experimental conditions presented before (Figures 1.2 (a) 

-(d)).: 


Sy =Zn(R' J /R j )//n(P , i/P i ), 


(7) 



= 0.01, Xo 2 = 0.001, Xat = 0.989, P av = 3.4 atm, ♦ for Xh 2 = 0.1, Xq 2 = 0.005, Xat = 0.895, P 
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where Pi={&i) and e(OH), and Rj = A max , NS max and tso- Figures 1.4 (a) - (d) 
show the sensitivities for A max , NS max and tso by raising {&J and e(OH) by 
30%. At 1234 K (Figure 1.4 (a)), reactions (1), (2), (3) and e(OH) are sensitive 
to NS max while reactions (1) and (2) are sensitive to tso- When the rate 
coefficients of reactions (1) and (2) were increased, the values of A^* were 
also increased because of more OH production, whereas the increase of ks 
caused the decrease of A max due to the consumption of OH by reaction (3), OH 
+ H2 = H2O + H. The sensitivity of P^nax on e(OH) is not shown here for obvious 
reasons. At 1556 K (Figure 1.4 (b)), the sensitivities for NS max and tso are 
basically the same as at 1234 K, while only reaction (1) shows a little 
sensitivity for A max . At Figure 1.4 (c) condition (1700 K, 3.2 atm), the 
sensitivities of NS max and tso ° n reaction (2) disappear and reaction (11) 
becomes important. Again reactions (1) and (3) are sensitive to A majc as in the 
case of 1234 K. At 2163 K (Figure 1.4 (d)), nothing shows sensitivity on A max ■ 
NS ma * is sensitive to reaction (1) and e(OH) while tso is sensitive to the 
reactions (1), (3) and (11). Overall, NS max is sensitive to reaction (1) and 
e(OH) while tso is sensitive only to reaction (1). Therefore, k\ could be 
determined by matching only characteristic times (t25, tso and t7s) without any 
interference from absorption characteristics. However, the effect of possible 
contaminants present must be considered. The influence of contaminants 
was investigated for test gas mixtures presumed to contain C3H8 or H-atom at 
the level of 0.05 to 0.5 ppm and 5 to 50 ppb, respectively. A reaction 
mechanism reported by Frenklach and Bornside [45] capable of simulating 
reaction paths for CH4 and CH4 / C3H8 mixture oxidation was used. Absolute 
sensitivities of tso defined as (Y'-Yo)/Yo, are given in Figures 1.5 (a)-(d). At low 
to middle temperatures (Figures 1.5 (a)-(c)), as Xc 3 h 8 increases tso decreases 
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Figure 1.4(a) Sensitivity spectrum for the experimental condition of Xh 2 - 
0.02, Xo 2 = 0.005, Xat = 0.975, T 5 = 1234 K and P 5 = 0.945 atm. Ordinate 
is logarithmic sensitivity defined in the text. Reaction numbers are listed 
in Table 1.3. Sensitivity was computed setting rate coefficients to 1.3 x 
Table 1.3 values and the absorption coefficient obtained in Figure 1.3. ■ 
NS,^* sensitivity, ED sensitivity and m ts 0 sensitivity. All noticeable 
changes are included. 
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Figure 1.4(b) Sensitivity spectrum for the experimental condition of Xh 2 = 
0.02, Xo 2 = 0.005, Xat = 0.975, T 5 = 1556 K and P 5 = 0.751 atm. Ordinate 
is logarithmic sensitivity defined in the text. Reaction numbers are listed 
in Table 1.3. Sensitivity was computed setting rate coefficients to 1.3 x 
Table 1.3 values and the absorption coefficient obtained in Figure 1.3. ■ 
NS^r sensitivity, E3 sensitivity and E tso sensitivity. All noticeable 
changes are included. 
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Figure 1.4(c) Sensitivity spectrum for the experimental condition of Xh 2 = 
0.01, Xo 2 = 0.001, Xa r = 0.989, T 5 = 1700 K and P 5 = 3.203 atm. Ordinate 
is logarithmic sensitivity defined in the text. Reaction numbers are listed 
in Table 1.3. Sensitivity was computed setting rate coefficients to 1.3 x 
Table 1.3 values and the absorption coefficient obtained in Figure 1.3. ■ 
NS,^ sensitivity, E3 P^nax sensitivity and E tso sensitivity. All noticeable 
changes are included. 
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Figure 1.4(d) Sensitivity spectrum for the experimental condition of Xh 2 - 
0.02, Xo 2 = 0.002, Xat = 0.978, T 5 = 2163 K and P 5 = 1.994 atm. Ordinate 
is logari thmi c sensitivity defined in the text. Reaction numbers are listed 
in Table 1.3. Sensitivity was computed setting rate coefficients to 1.3 x 
Table 1.3 values and the absorption coefficient obtained in Figure 1.3. ■ 
NS^uuc sensitivity, E3 A max sensitivity and Q tso sensitivity. All noticeable 
changes are included. 
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Figure 1.5(a) Sensitivity spectrum of tso by H-atom or C3H8 impurity for 
the experimental condition of Xh 2 = 0.02, Xo 2 = 0.005, Xat = 0.975, T 5 = 
1234 K and P5 = 0.945 atm. Ordinate is absolute sensitivity defined as (Y 
- Y 0 ) / Y 0 where Y 0 is the value of tso without H-atom or C3H8 impurity. E 
H-atom sensitivity and ■ C3H8 sensitivity. 
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Figure 1.5(b) Sensitivity spectrum of tso by H-atom or C3H8 impurity for 
the experimental condition of Xh 2 = 0.02, Xo 2 = 0.005, Xat — 0.975, T 5 = 
1556 K and P5 = 0.751 atm. Ordinate is absolute sensitivity defined as (Y' 
- Y 0 ) / Y 0 where Y 0 is the value of tso without H-atom or C3H8 impurity. E2 
H-atom sensitivity and ■ C 3 H 8 sensitivity. 
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Figure 1.5(c) Sensitivity spectrum of tso by H-atom or C3H8 impurity for 
the experimental condition of Xh 2 = 0.01, Xo 2 = 0.001 , Xat = 0.989, T 5 = 
1700 K and P5 = 3.203 atm. Ordinate is absolute sensitivity defined as (Y 
- Y 0 ) / Y 0 where Y 0 is the value of tso without H-atom or C3H8 impurity. ZZA 
H-atom sensitivity and ■ C3H8 sensitivity. 
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Figure 1.5(d) Sensitivity spectrum of tso by H-atom or C3H8 impurity for 
the experimental condition of Xh 2 = 0.02, Xo 2 = 0.002, Xa r = 0.978, T5 = 
2163 K and P5 = 1.994 atm. Ordinate is absolute sensitivity defined as (Y 
- Y 0 ) / Y 0 where Y 0 is the value of tso without H-atom or C3H8 impurity. EZ3 
H-atom sensitivity and ■ C3H8 sensitivity. 
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markedly, obviously owing to the increase of the C 3 H 8 decomposition rate to 
produce H-atoms. At 2163 K, H-atoms are generated by the fast initiation 
reactions (H + H + M = H 2 + M and H + HO 2 = H 2 + 0 2 ) so that contaminant 
effects were found to be small. At 1234 K, even though the effect of C 3 H 8 is 
not as big as at higher temperatures (1556 K and 1700 K) because of the slow 
decomposition of C 3 H 8 , the effect of H-atoms which could be produced easily 
from other possible contaminants (e.g., long-chain hydrocarbons such as 
vacuum pump oil), is great as seen in Figure 1.5 (a). As an example, matching 
tso, shortened about 10% by contaminants, requires about a 15% increase of 
k\. At all conditions tested the sensitivity of NS max to the added contaminants 
was negligible. 

Based upon the sensitivity analysis for {&il, e(OH) and contaminants, 
the parameters k\ and e(OH) were chosen for simultaneous optimization. 
Goodness of match to the experimentally determined NS ma * and A max was the 
criteria used in optimization. In the simulation, NS™^. the normalized slope at 
maximum OH growth, was generated from the time derivative of the light 
transmittance: 

NS max = ahs{d(I/Io) / dt } max 

= abs{-e(OH) x l x exp[-e(OH) x l x CohI * (dCoH / dt )} majc (8) 

where HI 0 is the transmittance, e(OH) is the absorption coefficient of OH, l is 
the light path length and Coh is the OH concentration. The NS mar is not only 
a function of e(OH) as seen in the above equation, but also being exclusively 
dependent upon k\. Thi s indicates that k\ and e(OH) are coupled together on 
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NS maa: and A max • This coupling can also be observed in the sensitivity 
analyses. For instance, the sensitivities of e(OH) on NS ma * amount to about 
half of the sensitivities of k\ at the conditions given before. Thus it seems 
difficult to separate the effects of k\ and e(OH) on NS ma * and A max unless a 
special statistical method is employed. However, if k\ is changed by 30% with 
a fixed value of e(OH), the sensitivities of NSma* are more than three times as 
large as those of A max . Therefore it is possible to simultaneously optimize k\ 
and e(OH) without relying upon the extensive parameter optimization 
techniques used elsewhere [9, 33, 46, 47, 48, 49]. In the simulation the 
absorption coefficient of OH was varied first to match A max , then k\ was varied 
to match NS ma3C • This procedure was repeated several times until 
satisfactory matches were achieved. The reaction rate coefficient expressions 
of k 2 and &3 in Table 1.2 were assumed to be correct and were not subjected to 
the optimization, although they showed some sensitivities at T < 1700 K. For 
k 2 , the two parameter rate coefficient expression of Sutherland et al. [31] (880 
K < T < 2495 K) was used. Recent results from this laboratory indicate that k 2 
could be best represented by an Arrhenius rather than a non-Arrhenius type 
expression in the temperature range of interest in this study. For k$, a recent 
fit to the existing data was taken [25] . 

The values of k\ obtained in this way are listed in the last column of 
Table 1. 1 and are presented in Figure 1.6. The least squares fit to the data for 
1050 K < T < 2500 K is given by 


) = (7.13 ± 0.31) x 10 13 exp(- 6957 ± 30 K/T) cm 3 mol' 1 s' 1 (9) 
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with uncertainties of 1 -ct standard deviation. In the simulation there was a 
temperature rise during the course of reaction. As k\ was determined using the 
OH growth early in the reaction and e(OH) was determined using the m a ximum 
absorption rate in the reaction, the temperatures at those relative reaction 
times varied by less than 25 K depending upon initial conditions. The 
expression for the temperature dependent shape function (g(co 0 )) for e(OH) 
given by Hanson and coworkers [50] gave negligible change in e(OH) during the 
reaction due to the increase of temperature. 

As the coupling between k\ and e(OH) is not reflected on {t25> tso and t 75 } 
and A max , k\ and e(OH) were also optimized using characteristic times and 
A m ax by ignoring possible contamination effect. The results are plotted in 
Figure 1.7. The solid line represents the least squares fit to the data and is 
given by 

fcidtil) = (7.19 ± 0.41) x 10 13 exp(- 7015 ± 40 K/T) cm 3 mol" 1 s" 1 (10) 

for 1050 K < T < 2500 K. The uncertainties are 1 -c standard deviation. The 1- 
c standard deviation of individual points from the ^i(NS ma3: ) and &i({ti}) are ± 
4% and ± 6% respectively. 

C. Error Analysis 

Error analyses for individual &i(NS max ) and &i({til)) were performed 
under high temperature conditions since the most probable uncertainties in the 
measurements of the incident shock velocity and of the pressure of the 
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Figure 1.6 Arrhenius plot of the experimental data of'^i(NS m a*)- The solid 
line is the least squares fit to the data, ^(NS,^*) = 7.13 x 10 13 exp(- 6957 
K/T) cm 3 mol" 1 s" 1 (1050 K < T < 2500 K). Symbols are A for 0 =2 
mix tures, □ for 0 = 5 mixtures, O for 0 = 10 mixtures, and • for 0 = 15 
mix tures (15.0% H 2 , 0.5% 02, 84.5% Ar). The data for 0 = 15 mixtures 
were not included in the least square fit for ^i(NS max ) (see text). 
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Figure 1.7 Arrhenius plot of the experimental data of &i({ti}). The solid 
hne is the least squares fit to the data, &i({ti}) = 7.19 x 10 13 exp(- 7015 
K/T) cm 3 mol _1 s _1 (1050 K < T < 2500 K). Symbols are A for =2 
mixtures, □ for <J> = 5 mixtures, O for <f> = 10 mixtures. 
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reflected shock region (P5) were from the high temperature experiments. The 
sampling time of the oscilloscopes used in the incident shock velocity 
measurements was set to 0.5 ps per point. If we consider 1 ps as the 
maximum uncertainty in time measurement, this uncertainty propagates to 
0.5% error in shock velocity, which yields 0.9% error in T5 and 0.5% error in k\. 
As described before, Ps ( used for the correction of T5, was obtained from the 
last pressure transducer signal. About 0.7% reading error of the signal output 
resulted in a 0.7% error for P 5 , which produced a 0.2% error in the corrected T 5 
and a 1.0% error in k\. The error in k\ from the initial test gas mixture 
pressure (Pi) was negligible. Although the possible errors introduced in the 
mixture preparation is small, the effect on is not ignorable and they are 
appropriately considered in the error propagation analysis below . 

The largest error introduced to &i(NS mo * ) was from NS max . The 
maximum error in the measurement of NS^^uc was 6%, which alone generated 
an error of 3% in /fei(NS max ). The matching error of less than 4% for the 
computed to the experimental NS max is within the error bounds of NS max 
considered above. Since the coupling between &i(NS max ) and A max through 
e(OH) is not discernible at high temperatures as shown in the sensitivity 
analysis, the error to ^lCNSma* ) from A max evaluation was examined under 
the low temperature condition ( 1234 K) of Figure 1.4 (a). Simulations indicated 
that a 3.0% error in A max could only change the value of &i(NS majc ) by less 
than 0.5%. 

The maximum uncertainty in taking characteristic times from 
experimental records was 2.0%, resulting in 2.0% error in &i({ti}). Because the 
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characteristic times are not coupled with e(OH), errors in A max do not induce 
any error in &i({til). The maximum error in matching, for example, tso, is 2.0% 
and the effect on &i({til) is also 2.0%. 

The propagation of error was computed using the following equation at a 
high temperature condition of Figure 1.4 (d), [51], 

s (F) = sqrtjf £w*l* + gfW + (f + • • • j (11 > 

where F is the individual data points of fei(NS m<T J or of &i({ti}) and x, y, z , ... are 
error sources. A very conservative error propagation analysis with all the 
possible errors mentioned above gives errors for ^i(NS max ) = ± 6.0% and for 
&i({til) = ± 7.5%. The percentage contributions to the errors for the data of 
&i(NS ma * ) are 54% for NS,^*, 24% for P5 in the temperature correction, 10% 
for incident shock velocity, and 10% and 2% for AXo 2 and AXh 2 of the initial 
mixture composition, respectively. A similar distribution of the error 
contributions was also obtained for £i({til). 


Discussion 

Figures 1.2 (a) - (d) show the computed profiles (smooth lines) by using 
the reaction mechanism in Table 1.2. At temperatures of 1556 K, 1700 K and 
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2163 K with different final pressures and mixture compositions, the Table 1.2 
reaction mechanism reproduced the experimental profiles quite well. However, 
as seen in Figure 1.2(a) (1234 K), the simulation underpredicts the OH 
concentration after the peak. Only more than a factor of 10 reduction of k% (H 
+ H + M = H 2 + M), which does not affect the pre-peak absorption profile at all, 
could bring the simulated OH concentration to the experimental level. This 
disagreement is possibly due to either the additional heating of the still moving 
gas toward the reaction zone long after the beginning of the reaction (~lms) so 
that the OH producing reactions proceeded faster while the computer 
simulations could not properly account for these phenomena, or it is due to the 
current incomplete understanding of HO 2 chemistry at these conditions. 
Despite this mismatch after the absorption peak at low temperatures, the k\ 
values reported in this study are valid because the rapid growth of OH 
concentration is mostly influenced by the chain branching reaction, H + O 2 = 
OH + O, and in any case the matches to the pre-peak experimental profiles are 
excellent. This is further manifested by the exact reproduction of the timewise 
change of the pre-peak OH concentration. 

The absorption coefficients of OH (Figure 1.3) show two distinct 
characteristics, i. e., dependence upon temperature and pressure. They show a 
small positive temperature dependence, absorption coefficients being slightly 
high at high temperatures. The temperature dependence comes mainly from 
the Boltzmann fraction in the integrated absorption (S) expression and the 
Voigt function (g(a> 0 ))- As temperature increases, the Boltzmann fraction 
decreases and the value of the Voigt function at fixed pressure increases 
because of the decrease of the spectral line broadening [27, 28, 29, 30]. 
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However, the positive temperature dependence indicates that the Voigt 
function is a stronger function of temperature than the Boltzmann fraction. 
The pressure dependence of e(OH) can be interpreted in terms of the spectral 
line broadening due to collision (collision broadening). The collision broadening 
width is directly proportional to the system pressure [27, 50] . As expected, 
e (OH) decreased as pressure increases at fixed temperature. Because the 
Voigt function is the convolution integral of Gaussian (Doppler) and Lorentzian 
(collision) function [27, 28, 29, 30], its value optimized for each experiment 
reflects the extent of contribution by collision broadening to e(OH) [52], 

In Figures 1.8 and 1.9 the matches between the experimental and 
computed NS max and A max are seen, within the data scatter of ± 4% and ± 
3%, respectively, to be excellent over the entire range of temperature and 
density of the various compositions of the mixtures studied. The mean 
deviations are 0.0 ± 0.96% for NS ma * and 0.15 ± 0.68% for A max . The 
uncertainties are 1-c deviation. The optimized A max could not be distinguished 
in the evaluation of £i(NS max ) and &i({til) in most cases. Figure 1. 10 shows 
the match of a characteristic time, tso. The mean value of the deviation is - 
0.42% with a standard deviation of 1.42%. The negative mean deviation of tso 
denotes that the simulated tso's are slightly shorter. However, it is within noise 
level, especially at high temperatures. For example, at 2163 K, the measured 
tso was 64.0 (is with an error of ± 1.0 |is or ± 2%. 

The &i({tj}) expression gives always slightly lower values than ^i(NS max ) 
over the full temperature ranges, but the differences are quite small; at 1050 
K, &i(NS ma * ) = 9.45 x 10 10 cnrhnol-V 1 , *i({ti}) = 9.02 x 10 10 cnr’mol-V 1 while 
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Figure 1.8 Match between experimental and calculated values of NS ma *. 
The percent deviation is defined as 100 x [NS max (calculation) — NS maa: 
(experiment)]/ NS^* (experiment). Symbols are A for <J> = 2 mixtures, □ 
for <|> = 5 mixtures, and O for <(> = 10 mixtures. 
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Figure 1.9 Match between experimental and calculated values of A^*. 
The percent deviation is defined as 100 x [h max (calculation) A// mx 

(experiment)]/ Pimax (experiment). Symbols are A for cf> = 2 mixtures, □ for 
0 = 5 mixtures, and O for 0 = 10 mixtures. 
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Figure 1. 10 Match between experimental and calculated values of tso. The 
percent deviation is defined as 100 x [tso (calculation) - tso (experiment)]/ 
tso (experiment). Symbols are A for 0 = 2 mixtures, □ for <J> = 5 mixtures, 
and O for 4> = 10 mixtures. 
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at 2500 K, ifeiCNS^ ) = 4.41 x 10 12 cm^oHs" 1 , *i({ti}) = 4.35 x 10 12 cmSmol" 1 
s _1 . If there were some impurities in the experimental system, the 
characteristic times would be shortened so that a higher &i({ti}) would be 
resulted in, which indicates that our system was not contaminated with 
impurities. As the above comparison indicates, ^i(NS maa; ) and ^i({tj}) overlap 
each other within the mutual error limits. This confirms that our experimental 
results were not influenced by any possible impurities present. 

The present determination of k\, Eq.(9), is compared to the previous 
experimental studies in Figure 1. 11. Yang et al. [13] performed a single beam 
laser absorption study of OH radical. A reinvestigation that surplanted their 
previous work on the title reaction [14, 15]. An iterative optimization using 
eight time - difference responses (e.g., tso-t^) and an adjustable time - zero 
yielded a rate coefficient expression that is 5% lower at 1850 K and 5% higher 
at 2500 K than our expression. Yang et al. then combined their data with that 
of Shin and Michael and obtained a new expression that is well within our error 
bounds, being 4% lower at 1100 K and 2% higher at 2500 K. 

Du and Hessler [11] measured the rate coefficients for reaction (1) in the 
temperature range of 3450 K to 5300 K by a flash absorption - shock tube 
technique. They combined their results with the data of Masten et al. [8], Shin 
and Michael [10] and Pirraglia et al. [7] and reported for the temperature range 
of 960- 5300 K: 


= (9.76 ± 0.72) x 10 13 exp(- 7474 ± 122 K/T) cnAnol-V 1 


(12) 
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Figure 1. 11 Comparison of the present results for ^(NS,^) to the previous 
experimental studies. 
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with uncertainties at the 95% confidence level. Comparing equations (9) and 
(12), the Du and Hessler's mean values are 16% lower at 1050 K and 11% 
higher at 2500 K than the present results. Since large error limits are involved 
in the above expression, the two expressions are in agreement within the 
combined experimental uncertainties over 1050K < T < 2500K. 

In their evaluation Du and Hessler introduced a profile shift parameter, 
x c , that was needed to shift their calculated profiles to longer distances in order 
to match the experimental profiles. They associated x c with the incubation 
time of H 2 dissociation, the main initiation channel at high temperature. 
However, we are able to model their Figure 3 experiment using the Table 1.2 
mechanism without resorting to the use of profile shifting. This comparison is 
shown in Figure 1.12, where the solid line is our predicted absorption profile and 
the dotted lines are the upper and lower noise band of the Du and Hessler's 
profile. For this condition Du and Hessler used a 0.783 mm shift of their 
calculated profile. Our prediction is well within the noise band, although 
approaching the high end at long distances, where the profile is controlled by H 2 
decomposition. Indeed, at all times, the profile is more sensitive to reaction (6) 
(H + H + M = H 2 + M) than all other reactions including reaction (1) (H + O 2 = 
OH + O). This can be seen in Figure 1. 12 where the upper filled boxes and the 
lower open boxes represent the extent of profile shift for ± 20% variations in k\ 
and ke respectively. 

Shin and Michael studied the reaction (1) by using the laser photolysis - 
shock tube technique by monitoring H-atom depletion in the pseudo-first-order 
conditions over the temperature range of 1103 - 2055 K [10]. The Arrhenius 
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Figure 1. 12. Comparison of our computed profile to the experimental profile 
of Du and Hessler's Figure 3 condition [11]. Solid line is the absorption 
profile computed using the Table II reaction mechanism with x c = 0 and the 
dotted lines are upper and lower noise band of the Du and Hessler's profile. 
Early part of experimental signal is dominated by shock passage 
transients. Filled and open boxes represent the extent of profile change for 

± 20% variation of &i(H + 02= OH + O) and &6(H + H + M = H 2 + M) 
respectively. 
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expression for k\ is given by 

*1 = (6.93 ± 0.96) x 10 13 exp(- 6917 ± 193 K/T) cm 3 mol- 1 s" 1 ( 13) 

where uncertainties are at the 1-c level. Since they derived k\ from the 
measured slope of ln( absorbance )*, their results are independent of impurities, 
just as our determination of k\ from the maximum slope (NS maa: ) was 
insensitive to impurity effects. It is important to emphasize that in Shin and 
Michael's experiments, corrections of T5 and ps for boundary layer effects were 
made using the same basis [23, 53] as was done in the present study. As seen 
in Figure 1. 11, the k\ values of Shin and Michael and of the present study are in 
quantitative agreement over the temperature range of 1103 - 2055 K. 

Masten et al. [8] determined k\ using laser absorption spectroscopy to 
detect OH radicals formed by heating mix tures of H2 and O2 diluted in Ar (<[> = 
-2.5, and -5.0 ) behind the incident and reflected shock waves for the 
temperature range 1450 - 3370 K. The Ri(5) line of the A 2 E + *- X 2 n system of 
OH at 306.687 nm (vacuum) was utilized to monitor OH concentration. The 
reported k\ expression is given by 

k\ — (9.33 ± 0.40) x 10 13 exp(- 7448 ± 86 K/T) cm 3 mol" 1 s' 1 ( 14) 

At the extreme ends of the common temperature range this expression 
is 7.5% higher at 2500 K and 7.0% lower at 1450 K than the present study; 
however the two experiments agree with each other within the given mutual 
error limits, ± 10% by Masten et al. and ± 6% by this study. Even at the 
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extreme ends of both temperature ranges, 3370 K (Masten et al.) and 1050 K 
(the present study), the k\ values agree within the combined error limits. 
Simulations were performed using the reaction mechanisms of Masten et al. 
and of Table 1.2. Laminar flow and constant volume conditions were assumed 
for incident and reflected shock wave experiments, respectively. At their 
Figure 3 condition (1980 K), maximum OH concentration and the 
characteristic times t 25 , tso and t 75 were simulated. Using either the Masten 
et al.'s or our k\ expression resulted in identical values for [OH] max . 
Furthermore, our mechanism reproduced the experimental results of Masten 
et al. as well as their reaction mechanism. However, our k\ value is 2% 
smaller than that of Masten et al. at 1980 K. Substituting our k\ expression 
into the reaction mechanism of Masten et al. and adjusting k\\ (H + H0 2 = H 2 
+ 0 2 ) to match characteristic times, as was done by Masten et al., a 34% 
increase of k\\ was required. Next, as there is no [OHl^o* under Masten et al.'s 
Figure 5 condition (2898 K), the first visible appearance of plateau OH 
concentration and the characteristic times, t 2 5 , tso and t 75 defined above were 
taken as simulation parameters. Again, our mechanism reproduced Masten et 
al.'s experimental results as well as the reaction mechanism of Masten et al. 
At this temperature, our k\ value is about 10% smaller than that of Masten et 
al. As before, we substituted our k\ expression and adjusted &n to obtain a 
match. A factor of 2.5 increase of k\\ was required. At this temperature, 
thermal decomposition rate coefficient of H 2 , k&, is an important initiation 
source, indeed, the sensitivities of ke and k\\ are equal. However, NS max is not 
sensitive to either k\\ or k% at this condition. In essence, the difference 
between the k\ expressions derived by Masten et al. and in this study are 
caused by secondary chemical reactions, i.e., either k\ expression can be used 
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to simulate the characteristic times of both sets of experiments by carefully 
choosing the rates of the initiation reactions. However, the temporal behavior 
of the OH concentration measured in this study was well reproduced using the 
k\ value determined by matching NS max . Moreover, this determination was 
independent of any cross sensitivity upon initiation reactions. The effect of 
and k% on ^iCNS^o*) were also examined under our conditions used above and 
was found to be negligible. Table 1.2 reaction mechanism with and ks 
expressions of Masten et al. resulted in 0.6% and 1.3% changes in k\ values at 
1234 K and 2163 K though there are 3 - 20% and 3-6% differences in k<i and k% 
respectively. 

Yuan et al. published a k\ expression determined by laser absorption of 
OH radicals generated by heating mixtures of hydrogen and oxygen diluted in 
Ar behind reflected shock waves over the temperature range 1050-2700 K [9]. 
Five active parameters, k\, k% ks, k-n, and AfH°298(H02) were simultaneously 
optimized by minimizing the difference between the computed and 
experimental characteristic times, t 25 , tso, t 75 , using the solution mapping 
method [54] . The reported non-Arrhenius type of k\ expression is 

ki = 1.59 x 10 17 T~°- 927 exp(- 8493 K/T) cm 3 mol' 1 s” 1 ( 15) 

with oiogAj = 0.05. Their result is in excellent agreement with that of Schott [3] , 
1.22 x 10 17 'T 0 907 exp(- 8369 K/T) for 1250 K < T < 2500 K. Compared to the 
present experimental results in the common temperature range, 1050 K < T < 
2500 K, the k\ values of Yuan et al. are on average about 17% smaller, but 
agree with each other within the combined error limits. These authors also 
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reported that with Masten et al.'s secondary reaction rate coefficients, k 2 , h, 
and /fen, a slightly different expression was resulted, 1.71 x 10 17 T -0 - 932 
exp(-8498 K/T), which differs from the above equation (17) by -3%. The effect 
of the secondary reaction rate coefficients on the determination of k\ was thus 
not large (also see Ref. 14, 15). 

Recently the data of Yuan et al. and Masten et al. were subjected to a 
rigorous optimization using two different model responses by Yu et al. [12], the 
characteristic times (t 25 , tso, t 7 s) and the time difference. At, defined as At = t 75 
- 125- The data of Masten et al. yielded nearly identical k\ expressions for both 
model responses while those of Yuan et al. were noticeably different, with 
^(characteristic times) being 15% lower than &i(A() at 2150 K. There was 
reasonable agreement between the Yuan et al.'s ki(A t) and both the Masten et 
al.'s k\ expressions (which were nearly identical). An explanation for the 
difference between the Yuan et al.'s k\(characteristic times ) and &i(A£) was 
proposed based upon the effect of vibrational non-equilibrium of O 2 . The high 
temperature data of Yuan et al. was taken at short reaction times where this 
effect would be most pronounced. This is akin to the explanation given by Du 
and Hessler for the profile shift parameter included in their study. Yu et al. 
suggest an expression based upon the At response surface that is within our 
error bounds, being 3% higher at 2500 K and 7% lower at 1336 K than the 
present study. 

The possible effect of O 2 vibrational relaxation on hydrogen combustion 
has important implications for the performance of airbreathing hypersonic 
propulsion systems. Such systems would be hydrogen fueled and have 
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characteristic chemical reaction times much shorter than 100 ps. In order to 
quantify the effect of O 2 vibrational relaxation a series of experiments was run 
using the Yuan et al.'s Series E composition in the pressure and temperature 
range where the effect should be most visible, i.e. 2.4 to 3.0 atm and 1905 to 
2380 K, respectively. Unlike the case of the Yuan et al.'s data we did not see a 
difference between the k\ value determined using either tso or At as the 
modeling criteria. Experimental values of tso and A£ are shorter than predicted 
using the Table 1.2 mechanism, on average by 2% (0.6 (is) and 5% (0.3 ps) 
respectively. While arbitrary time accuracy may be achieved in the 
simulation, differences of this scale, although they are discernible, are not 
particularly meaningful. Nonetheless, we can match tso, and NS max for 
each experiment. At these conditions there is exquisite sensitivity to tso as 
modeling criteria - a 0.1 ps change required a 1% change in k\, while sensitivity 
to At was smaller. The optimized k\ values for these experiments differ on 
average 4% (within our uncertainty limits) from our recommended expression 
(individual experiments required adjustments between - 3% and 8%). We are 
not required to invoke O 2 vibrational relaxation to explain our short reaction 
time data as both tso and At are predicted equally well using our mechanism 
without additional constraints. It is worth remembering that Belles and 
Lauver [55] had previously shown that O 2 vibrational relaxation was not 
required to explain the lengthening of induction delay at short times in H 2 /O 2 
mixtures that had been proposed by Schott and Kinsey [56] . 

It is worthy to note the flash photolysis - shock tube work of Pirraglia et 
al. in the temperature range 962 to 1705 K, because it has been used in all of 
the recent studies either for comparison^, 9, 10, 15] or to make fits for k\ in 
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the form of non- Arrhenius type [ 8 , 15] or Arrhenius type [11] by combining 
with their own results. In the experiments, a small amount of either NH 3 or 
H 2 O was photolyzed in excess of O 2 and Ar and H-atom decay was monitored 
by ARAS. The k\ expression reported is 

k\ = (1.68 ± 0.19) x 10 14 exp(- 8119 ± 139 K/T) cm 3 mol -1 s _1 (16) 

with the mean deviation of ± 16%. The above expression for k\ is compared to 
our results. The two sets of expressions agree within ± 22% over the common 
temperature interval with k\ of the present study being 22% larger at 1050 K 
and 19% smaller at 1700 K. However, as pointed out by one of the authors in a 
recent review [32], the k\ expression of Pirraglia et al. was influenced by the 
data at low temperatures, T < 1200 K. Careful inspection of the data (Figure 2 
of Ref. 7) showed that there seemed to be two sets of data separated at 1200 K 
and that the above k\ expression overpredicts the data at T < 1200 K. These 
two sets of data could be separately fitted into 6.76 x 10 13 exp(- 6829 K/T) cm 3 
mol " 1 s _1 (1206 K < T < 1705 K) and 1.56 x 10 14 exp(- 8076 K/T) cm 3 mol ' 1 s ' 1 
(962 K < T < 1206 K). The former expression agrees with our results within 
better than 6 % in the common temperature range (Figure 1.13). The possible 
reason for the low k\ values at lower temperatures may be due to the over- 
correction either for temperatures for the boundary layer effect or for the 
contribution from reaction ( 9 ), H + O 2 + M= HO 2 + M. It is also possible that 
the k\ values actually start to decrease at temperatures below 1200 K so that 
the k\ expression has a negative A-factor. Nevertheless, the upward 
temperature corrections ranged from 4% to 6 % of the ideal shock 
temperatures (average 50 to 60 K) for experiments with initial pressures of 



% Deviation log( k x /cm 3 mol 


11.0 


lOkK/T 


Figure 1. 13 Comparison between our results for ^i(NS ma x) with the data of 
Pirraglia et al. [Ref. 7]. (a) Arrhenius plot, (b) The percent deviation defined 

as 100 x [k i( experiment) - &i(fit, this study)]/ &i(fit, this study). Symbols 
are O for this study and • for Pirraglia et al. 
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about 10.5 Torr and 15.5 Torr. The range of the temperature correction will 
vary depending upon shock tube characteristics and experimental conditions. 
However, a 50 to 60 K correction seems too large. If temperatures were over- 
corrected, then the actual temperatures would be lower so that the A-factor 
and the temperature dependence of the Arrhenius expression for k\ would 
decrease. As temperatures are lowered by reducing over-correction, the 
contribution from reaction (9) to the measured overall H-atom decay 
constants from which k\ values were obtained may not be ignorable. However, 
the influence of reaction (9) on the evaluation of k\ values is rather small 
compared to increasing T -1 Cowering the temperatures). 

The temperature dependence of the A-factor was reported by several 
authors: -0.907 by Schott, -0.927 by Yuan et al., -0.7 by Masten et al., 2.0 by 
Fujii et al. from the combination of their results with the data of Pirraglia et al.; 
and -0.816 by theoretical calculation of Miller [17, 18]. Miller [17] calculated 
the rate of the reaction (1) using quasi classical trajectory and quantum 
mechanical threshold (QCT-QMT) methods and the potential energy surface of 
Melius and Blint [19], and attributed the negative temperature dependence to 
the nonstatistical "recrossing" effects especially at high temperatures. 
However, the results of the present work and Shin and Michael do not support 
this negative temperature dependence for HOOK < T <2500K. Recently 
Varandas et al. [20] calculated the thermal rate coefficient of reaction (1) using 
the quasiclassical trajectory (QCT), quasiclassical trajectory-quantum 
mechanical threshold (QCT-QMT), quasiclassical trajectory-internal energy 
quantum mechanical threshold (QCT-IEQMT), and the improved version of 
QCT-IEQMT methods at five temperatures, 1000, 1750, 2000, 2500, and 3000 
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K. In the calculation, the fourth version of their own double many-body 
expansion (DMVE IV) potential energy surface for the ground state of HO 2 , 
which reproduces the most accurate estimates of the experimental 
dissociation energy, equihbrium geometry, and quadratic force constants [57], 
was utilized. The results from the QCT-QMT methods gave best agreement 
with the present study within computational error bounds and are plotted in 
Figure 1. 14. Our k\ value is 4% smaller at 3000 K and 12% smaller at 1000 K 
where the minimum number of trajectories was calculated because of the 
computational burden. The Arrhenius expression for k\ from QCT-QMT 
methods given by Varandas et al. is k\ = (7.19 ± 1.11) x 10 13 exp(- 6943 K/T), 
which is very similar to equation (9) of the present study. 

There exists reasonable agreement between this study and most of the 
recent evaluations of k\ [7, 10-13, 58] so it is possible to achieve a consensus 
expression, shown in Figure 15, that is given by: 

fci = 7.82 x 10 13 exp(- 7105 K/T) cm 3 mol" 1 s" 1 (17) 

over the temperature range 960 to 5330 K, with an uncertainty of 6%. We 
developed this expression in the following fashion. The expressions of Pirraglia 
et al., Shin and Michael, Du and Hessler, Yu et al., Yang et al. and the present 
study were converted to a series of "data points" evenly spaced in 1/T over 
each of their temperature ranges. A weighted least squares fit was then 
obtained with the weighting factor being the inverse of the uncertainty limits 
for the Arrhenius expressions reported in the individual studies. The data of 
Pirraglia et al. and Shin and Michael were used to develop the Du and Hessler's 



Lg 


lOkK/T 


Figure 1.14 Comparison of our results for ii(NS mox ) with the 
experimental results of Masten et al. [Ref. 8] and the theoretical 
calculation of Varandas et al. [Ref. 20]. The solid line represents 
&i(NS mo *) expression. Symbols are A for Masten et al. and • for 
Varandas et al. 
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Figure 15. Consensus expression for k\. Solid line is the consensus 
expression, fciCNS,^) = 7.26 x 10 13 exp(- 6985 K/T) cm3 mo l-i s -i (960 K 
< T < 5330 K). Symbols are: V, Pirraglia et al.; ffl, Du and Hessler; □, Shin 

and Michael; A, Yu et al.; O, Yang et al.; O, This Study. Uncertainty 
limits for the individual expressions are shown for the highest and lowest 
temperature points. 
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expression and the data of Shin and Michael were used to develop the Yang et 
al.'s expression. Accordingly, their weighting factors were reduced to correct for 
this over representation. The weighting factors used to extract a consensus 
expression were 2 for Shin and Michael, 14 for Yu et al., 3 for Pirraglia et al., 8 
for Yang et al., 5 for Du and Hessler, 40 for this study, respectively. An 
Arrhenius expression was assumed. 


Conclusions 

The rate coefficient of reaction (1), H + 0 2 = OH + O was determined by 
using the OH laser absorption spectroscopic method behind reflected shock 
waves over the temperature range 1050 - 2500 K and pressure range from 
about 0.7 atm to about 4.0 atm. Eight different compositions of mixtures of <t> 
= 2, 5 and 10 were used in the study. Two different methods utilizing the 
absorption maximum (A max) and the normalized ma ximum slope (NS^*) and 
Amax and the characteristic times ({til) employed in the evaluation of k\ yielded 
almost identical results. The recommended expression for k\ is from the 
normalized maximum slope analysis and is given by 

ki = (7.13 ± 0.31) x 10 13 exp(- 6957 ± 30K/T) cm 3 mol" 1 s" 1 (18) 

with uncertainties of 1 -ct level. This expression is in quantitative agreement 
with that of Shin and Michael and is very similar to the computational results 
of Varandas et al. A complete comparison was given between the results of 
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the present study and those of recent experimental studies. The experimental 
data for k\ were best fitted using the Arrhenius type rate coefficient 
expression. Therefore, in the temperature range of the current study, which is 
relatively wide, we support the Arrhenius type expression for k\ although more 
data, especially below 1000 K, are needed to confirm the curvature. 

Critical review of the recent k\ determinations yields a consensus 
expression given by: 

ki = 7.82 x 10 13 exp(- 7 105 K/T) cm 3 mol" 1 s" 1 ( 19) 

over the temperature range 960 to 5300 K. 

The inverse reactant density dependence of k\ at 2000 K, which was 
suggested by Schott, could not be found even though the reactant densities 
were varied ten-fold. 
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Chapter II. O + H 2 = OH + H Reaction 


Introduction 

The reaction between oxygen atoms and hydrogen molecules 

O + H 2 = OH + H ArlTzsg = 7.8 kJ/mol (2) 

which is slightly endothermic, is the second most important chain branching 
step after the reaction of 

H + 0 2 = OH + 0 A r H °298 = 70.2 kJ/mol (1) 

in the H 2 /0 2 system. Although it was reported that flame propagation and 
ignition were not sensitive to the rate coefficients of reaction (2) [1], reaction 
(2) fur ni shes a reactive chain center (H-atom) directly, and also indirectly to 
the primary chain branching step, reaction (1), through 

OH + H 2 = H 2 0 + H A r H °298 = - 62.8 kJ/mol . (3 ) 

Because of the reason stated above, reaction (2) has been an interesting 
subject of research for many authors. Numerous experimental and theoretical 
investigations have been performed to determine the rate coefficients of 
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reaction (2). It is generally well known that the rate coefficients of this 
reaction, k 2 , can be represented by a curved non- Arrhenius type expression. 
The reaction rates are enhanced due to quantum mechanical tunneling [2] at 
low temperatures and due to the participation of the vibrationally enhanced 
reaction, 0( 3 P) + H 2 (u > 1) [3] at high temperatures. However, the extent of 
the curvature, especially, the vibrational enhancement at high temperatures 
is still not known well. 

Currently there are three reviews available for the rate coefficient 
information of reaction (2). Low temperature data obtained using discharge- 
flow systems and a few high temperature shock tube data up to 1972 were 
reviewed by Baulch et al. [4], The recommended k 2 expression in the 
temperature range of 400 - 2000 K is 

k 2 = 1.8 x 10 10 T 1 0 exp(- 4480 K/T) cm 3 mol -1 s" 1 (4) 

with the suggested error limits of ± 30%. Cohen and Westberg [5] reviewed k 2 
data available up to 1983 and recommended 

k 2 = 1.1 x 10 4 T 2 8 exp(- 2980 K/T) cm 3 mol -1 s _1 (5) 

for temperatures from 298 to 2500 K. The uncertainty was <Jiog(fe 2 ) = — 0-3 
298 K and at 1600 K and above. Because of the uncertainty and scatter in the 
data above 900 K, they suggested that an Arrhenius expression, k 2 = 7 x 10 13 
exp(- 5400 K/T) cu^mol^s" 1 described the data better at high temperatures (T 
> 1000 K). Warnatz [1] recommended the following k 2 expression which 
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differed only slightly from that of Baulch et al. (300 K < T < 2500 K). 

k 2 = 1.5 x 10 7 T 2 0 exp(- 3800 K/T) cm 3 mol* 1 s" 1 (6) 

Recent experimental results and one theoretical estimation for k 2 are 
listed in Table II. 1. 

From the rate coefficient expressions in the three reviews and in Table 
II. 1, it is easily seen that the k 2 values at combustion temperatures are widely 
different. For example, at 2000 K they sire (in units of cm 3 mol' 1 s" 1 ): 8.97 x 
10 12 (Warnatz, Ref. 1); 3.83 x 10 12 (Baulch et al., Ref. 4); 4.34 x 10 12 (Cohen 
and Westberg, Ref. 5); 6.93 x 10 12 (Schott et al., Ref. 6); 7.24 x 10 12 
(Pamidimukkala and Skinner, Ref. 7); 6.07 x 10 12 (Sutherland et al., Ref. 8); 
6.80 x 10 12 (Sutherland et al., extended temperature fit, Ref. 8); 7.04 x 10 12 
(Natarajan and Roth, Ref. 9); 6.55 x 10 12 (Natarajan and Roth, extended 
temperature fit, Ref. 9); 7.25 x 10 12 (Shin et al., Ref. 10); 6.89 x 10 12 (Davidson 
and Hanson, Ref. 11); and 7.74 x 10 12 (Garrett et al.. Ref. 12). The 
discrepancies worsen as the temperature increases. Also, in the Arrhenius 
expressions at high temperatures, the values of the A-factor (cm 3 mol _1 s _1 ) 
range from 1.87 x 10 14 to 8.13 x 10 14 and the temperature dependence from 
6854 to 9540 K, respectively. Furthermore, in the extensive study by 
Sutherland et al. [8], the difference between the values from the high 
temperature fit and the extended temperature fit, which incorporated low 
temperature data, becomes large as temperature increases. This invokes the 
question about the extent of curvature of the non-Arrhenius expression of k 2 . 



79 


Table II. 1. Results of the previous studies for k 2 a 


A (cm 3 mol *s 

n 

6(K) 

T range (K) 

Reference 

2.20 x 10 14 

0.00 

6916 

1400 - 1900 

6 

2.30 x 10 14 

0.00 

6916 

1919 - 2781 

7 

1.87 x 10 14 

0.00 

6854 

880 - 2495 

8 

4.34 x 10 13 

0.00 

5249 

504- 923 

8 

5.06 x 10 4 

2.67 

3165 

297 - 2495 

8 b 

3.72 x 10 6 

2.17 

4080 

1713 - 3532 

9 

3.87 x 10 4 

2.70 

3150 

300 - 3530 

9° 

7.90 x 10 14 

0.00 

9382 

1790 - 2250 

10 

8.13 x 10 14 

0.00 

9540 

2120 - 2750 

11 

2.89 x 10 14 

0.00 

7240 

1400 - 1900 

12 


a Rate coefficient expression is in the form of k 2 = A T n exp(- 0 / T). 

b Extended temperature fit using the results of their own and low temperature 
data (277 - 471 K) of Presser and Gordon [13]. 

c Extended temperature fit using the results of their own and those of 
Sutherland et al. [8] and low temperature data (277 - 471 K) of Presser 
and Gordon [13]. 
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In this study a series of experiments was performed at high 
temperatures in order to answer the two questions described above, i.e., what 
the absolute values of ^2 are; and how adequate the non- Arrhenius type 
expression for k<i is. 


Experimental Section 

The experimental set up and methods are the same as described in 
Chapter I. A mix ture composition of Xh 2 =0.0005, Xo 2 =0.005, and Xat=0.9945 

was selected for investigation based on the sensitivity analysis (see Computer 
Simulation). 


Results 

A typical hydroxyl radical absorption profile is shown in Figure II. 1. The 
rapid rise of the signal after an induction period is due to the exponential growth 
of OH by the chain branching and chain propagating reactions. Depending 
upon experimental conditions the OH concentration reaches a partial 
equilibrium or a maximum. 

As in the previous study of the H2/02/Ar system for k\ (see Chapter I), a 
set of experimental observables, {tso, P^max and NS ma *} were obtained from 
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Figure II. 1 A typical experimental record of OH absorption for the 
experimental condition of Xh 2 = 0.0005, Xo 2 = 0.005, Xat = 0.9945, T 5 = 
1898 K and P 5 = 1.826 atm. Characteristic time (tso), A max , and NS max 
were chosen as simulation parameters. The initial peak is the schlieren 
signal corresponds to passage of the reflected shock front. The smooth line 
indicates the computed OH absorption profile using the Table II.3 reaction 
mechanism and the OH absorption coefficient obtained from self- 
calibration at the maximum absorption (-Ana*). 
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each experimental record to represent the rate of the reaction progress and 
used for optimization targets of k^- The definitions of the observables are: A maa: 
is the absorption maximum, (1 - I/I o)max, tso is the characteristic time at which 
the light absorption has reached 50% of A max , NS max is the normalized 
maximum slope of OH growth, respectively. 

Table 1 1. 2 shows the experimental conditions and results. The 
temperature ranged from 1424 K to 2427 K and the pressures were kept 
around 1.9 atm. Temperature corrections for the boundary layer effects 
always increased T 5 from the ideal shock temperature by an average about 1.4 
%. Contrary to the results reported by Michael and Sutherland [14], the 
dependence of the ideal shock temperature and the initial pressure (Pi) upon 
the extent of T 5 correction could not be observed. 

A. Computer Simulation 

Computer simulations were performed using the detailed reaction 
mechanism and rate coefficient expressions in Table II. 3. The reverse reaction 
rate coefficients were automatically computed from the principle of detailed 
balancing. NASA thermodynamic data [15] were used to calculate the shock 
parameters and the equilibrium constants. The LSODE integrator [16] was 
utilized to solve a set of stiff differential equations describing chemical kinetics 
under assumed constant density conditions, calculated for the reflected shock 
waves [17], 
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Table II. 2 Experimental Conditions and Results of k<f 


t 5 

Ps 

Amax 

NS^ajc 

t50 

k 2 b ! 10 12 

fc 2 c / 10 12 


Xh 2 = 

=0.0005, Xo 2 = 

=0.005, Xat=0.9945 (<>=0.05) 


1424 

1.801 

0.198 

850 

797 

1.50 

1.50 

1430 

1.813 

0.190 

865 

795 

1.55 

1.58 

1449 

1.837 

0.198 

951 

756 

1.54 

1.62 

1528 

1.903 

0.238 

1444 

558 

2.28 

2.14 

1529 

1.912 

0.238 

1448 

592 

1.96 

1.99 

1559 

1.967 

0.235 

1579 

505 

2.29 

2.36 

1564 

1.964 

0.263 

1725 

508 

2.23 

2.30 

1671 

1.829 

0.259 

1990 

422 

2.94 

3.04 

1691 

1.844 

0.247 

2086 

382 

3.25 

3.31 

1715 

1.887 

0.220 

1774 

371 

3.19 

3.19 

1732 

1.942 

0.257 

2258 

335 

3.58 

3.40 

1769 

1.898 

0.295 

2601 

322 

3.65 

3.71 

1898 

1.826 

0.269 

2685 

252 

5.00 

4.70 

1948 

1.903 

0.280 

3260 

225 

5.26 

5.59 

1951 

1.900 

0.259 

2925 

219 

5.29 

5.40 

2093 

2.008 

0.294 

4215 

160 

7.14 

7.36 

2123 

1.928 

0.305 

4040 

157 

7.40 

6.79 

2127 

1.935 

0.306 

4271 

157 

7.44 

7.77 

2285 

2.004 

0.318 

5078 

125 

9.12 

8.83 

2334 

1.950 

0.299 

4886 

118 

9.90 

9.90 

2350 

1.964 

0.300 

5018 

114 

10.1 

9.94 

2370 

1.993 

0.305 

5374 

109 

10.4 

10.8 

2405 

2.050 

0.297 

5339 

101 

10.9 

10.9 

2416 

2.041 

0.301 

5281 

98 

10.8 

10.3 

2427 

2.064 

0.296 

5382 

94 

11.1 

11.3 


a Units are K for T 5 , atm for P 5 , |is for iso and cm 3 mol 1 s 1 for k 2 - 

b k 2 (t$o) 
c ^(NS,^*) 
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Table II. 3 Reaction Mechanism for k<f 



Reaction 

A 

n 

0 

Reference 

1. 

H + 0 2 = OH + 0 

7.13E+13 

0.0 

6957 

Chapter I 

2. 

0 + H 2 = OH + H 

1.88E+14 

0.0 

6897 

This study 

3. 

OH + H 2 = H 2 0 + H 

2.14E+08 

1.52 

1736 

19 

4. 

0 + H 2 0 = OH + OH 

4.51E+04 

2.70 

7323 

20 

5. 

o + o + m = o 2 + m 

Ar=1.0, H 2 =2.9, 0 2 =1.2, H 2 0= 

1.00E+17 

=18.5 

-1.0 

0 

21 

6. 

h + h + m = h 2 + m 

Ar=1.0, H 2 =4.0, H 2 O=12.0, H 

6.40E+17 

=26.0 

-1.0 

0 

21 

7. 

H + 0 + M = OH + M 
Ar=1.0, H 2 O=5.0 

6.20E+16 

-0.6 

0 

22 

8. 

H + OH + M = H 2 0 + M 
Ar=1.0, H 2 =2.5, H 2 0=16.25 

8.40E+21 

-2.0 

0 

21 

9. 

H + 0 2 +M = H0 2 + M 7.00E+17 

Ar=1.0, H 2 =3.33, 0 2 =1.33, H 2 0=21.3 

-0.8 

0 

21 

10. 

ho 2 + H = OH + OH 

2.20E+14 

0.0 

710 

23 

11. 

ho 2 + h = h 2 + 0 2 

2.50E+13 

0.0 

350 

1 

12. 

ho 2 + H = h 2 o + 0 

5.00E+12 

0.0 

710 

23 

13. 

H0 2 + 0 = 0 2 + OH 

2.00E+13 

0.0 

0 

1 

14. 

H0 2 + OH = H 2 0 + 0 2 

2.00E+13 

0.0 

0 

1 

15. 

HO 2 + HO 2 = H 2 O 2 + O 2 

1.06E+11 

0.0 

-855 

24 

16. 

h 2 o 2 + m = oh + oh + m 

Ar=0.67, O 2 =0.78, H 2 O=6.0 

1.20E+17 

0.0 

22900 

21 

17. 

h 2 o 2 + h = ho 2 + h 2 

1.70E+12 

0.0 

1900 

1 

18. 

H 2 0 2 + H = H 2 0 + OH 

1.00E+13 

0.0 

1805 

1 

19. 

H 2 0 2 + 0 = H0 2 + OH 

2.80E+13 

0.0 

3225 

24 

20. 

H 2 0 2 + OH = h 2 o + ho 2 

7.00E+12 

0.0 

720 

1 


a Rate coefficients are in the form & = A T n exp(- 0 / T). 
Units are cm 3 , K, mol, and s. 
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The local logarithmic response sensitivities defined below [18] was 
computed at a typical mid-temperature experimental condition presented in 
Figure II. 1. 


Si j =Zn(R , j /Rj)//n(P , i /Pi) (7) 

where P i ={>fe 1 } and R, = {tso and NSJj. Figure II.2 shows the sensitivities for 
t 5 o and NS ma x by raising the by a factor of 2. Reactions (1), (2), (3), and 
(11) were sensitive to characteristic time and reactions (1), (2), (3) and (4) 
were sensitive to NS max . Similar sensitivities were also observed for other 
reaction conditions. Even if reaction (2) has the highest sensitivity, it was 
difficult to isolate the effect of reaction (2) only on the changes of 
characteristic time and NS^. For the very lean condition (for example, $ = 
0.025), isolation of the contribution from reaction (2) could be obtained 
because the sensitivity of reaction (2) for tso and NSma* was about twice as 
big as those of the next highest one. However, in this lean mixture tso was 
extremely prone to contaminant effects. Therefore, the Xh 2 =0.0005, Xo 2 = 
0.005, X Ar =0.9945 mixture (<|> = 0.05) was chosen to study k 2 . In the 
computer simulation, k 2 values and e(OH) were varied to match tso, ^max and 
NS^o* ass uming that expressions of k\ from our previous extensive study, ks 
and &4 from the recent critical reviews by Oldenborg et al. [19] and Michael 
[20], respectively, and &n recommended by Warnatz [1] were correct. 
Especially, the effect of k n on the evaluation of k 2 was thoroughly searched 
because there are large discrepancies among the existing data at high 
temperatures (see Discussion). As was pointed out in Chapter I, the coupling 
between [til and e(OH) was found to be small. Therefore, first, A max was 
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0.5 


0.0 


- 0.5 


Figure II.2 Sensitivity spectrum for the experimental condition of Xh 2 = 
0.0005, Xo 2 = 0.005, Xat = 0.9945, T 5 = 1898 K and P 5 = 1.826 atm. 
Ordinate is logarithmic sensitivity defined in the text. Sensitivities are for 
200% increase to Table II. 3 values. Reaction numbers are listed in Table 
II.3. ■ for tso sensitivity and EZ3 for NS,^* sensitivity. 
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matched with e(OH) and then tso and NS max were matched by varying k 2 only. 

The values of k 2 obtained using tso and NS^x are listed in the last two 
col umns of Table II. 2 and are presented in Figure II. 3. The least squares fit to 
the combined data gives 

k 2 = (1.88 ± 0.07) xlO 14 exp(- 6897 ± 53 K/T) cm 3 mol -1 s -1 (8) 

with the uncertainties of 1- o standard deviation. The 1- c level deviation of the 
individual points from the fitted line is ± 6%. 

B. Error Analysis 

The sampling time of the oscilloscopes used in the shock speed 
measure m ent was set to 0.5 [is per point. Thus the maximum allowable 
uncertainty in the time measurement is 1.0 ps. An 1.0 [is uncertainty causes 
a 0.5% error in the shock speed, which yields a 0.9% error in T5 and a 0.5% 
error in k 2 . As described previously, P5, used in the correction of T5, was 
obtained from the last pressure transducer signal. A 0.7% reading error of the 
voltage output gave a 0.7% error for P5, which in turn produced a 0.2% error in 
the corrected T5 and an 1.0% error in k 2 . The error in k 2 from the m a xim u m of 
± 0.2% reading uncertainty of the initial pressure (Pi) was negligible. Although 
the possible errors introduced in the mixture preparation were small, the effect 
on k 2 was not ignorable. They were appropriately considered in the error 
propagation analysis. The errors in taking tso and NS ma x from the 
experimental records were 2.0% and 6.0%, respectively, and these errors in 



m 


13.0 



12.0 


lOkK/T 


Figure II. 3 Arrhenius plot of the experimental data of & 2 . O are 
Jfe 2 (NS max ), □ are ^ 2 ( 4 o). the solid line is the least squares fit to the 

combined NS max and £50 data, fc 2 = 1.88 x 10 14 exp(- 6897 K/T) cm 3 mol _ 1 s _1 

(1420 K<T< 2430 K). 
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turn produced 5.0% and 6.0% errors in k 2 , respectively. As seen in Figure II.4 
and Figure II.5, the m aximum errors in matching both tso and NSma* are less 
than 3%. 

Again, the error propagation was computed using equation ( 12) in the 
Chapter I. Considering all the errors, including the data scatter and the 
matching error, a maximum allowable error of ± 10% results for the k 2 

expression. 


Discussion 

A computed profile (smooth line) is shown in Figure II. 1. The reaction 
mechanism and the rate coefficient expressions in Table II. 3, especially the 
combination of k\, k 2 , k 3 , & 4 and k\\, reproduced the experimental profile quite 
well. The same level of matching is obtained for other experimental conditions. 
In Figures II.4 and II. 5 the matches between experimental and computed tso 
and NS ma * are shown, within the scatter of ± 3.0% for both t 5 o and NS ma *- 
Considering 2.0% and 6.0% for extracting tso and NS mar , respectively, from 
experimental records, the matches are excellent. In the optimization of k 2 
using tso and NS max , A max was always matched exactly to the experimental 
values. Matching A mnY resulted in a self-calibration of the absorption profile. 

The present determination of k 2 , Eq.(8), is compared to the previous 
studies in Figure II.6. Schott et al. [6] determined the ratio of k 2 /kx by 
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Figure II. 4 Matches of the computed tso to the experimental data. In 
computation, Table II. 3 reaction mechanism and rate coefficients with 
values in Table II. 3 were vised. The percentage deviation is defined as 100 

x {t 5 o(calculation) - ts 0 ( experiment)} / ^(experiment). 
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Figure II.5 Matches of the computed NS^ to the experimental data. In 
computation, Table II. 3 reaction mechanism and rate coefficients with 
values in Table II. 3 were used. The percentage deviation is defined as 100 

x {^^(calculation) - NS^( experiment)} / NS^Cexperiment). 
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Figure II. 6 Comparison of the present results for k<i to the previous 
studies. 
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measuring the O-atom concentration during the incident shock-wave initiated 
combustion of rich H 2 /O 2 /CO/CO 2 /AJ mixtures (<j> = 3.0, 5.0 and - 8.1) for the 
conditions of 1400 K < T 2 ^ 1900 K and 0.25 atm < P 2 ^ 1.73 atm. Light 
emission intensity from CO + O — » CO 2 * + hv was monitored. [O] was 
calibrated against the partial equilibrium transmitted light intensity, I pe and 
[OIpe, obtained with H 2 /O 2 = 1 mixture. The spike O-atom concentration, [0]s, 
calculated from each experiment with the aid of calibration experiments 
mentioned above, was simulated by varying kil k\. The results are presented 
as 


&2 = 2.2 x 10 14 exp(- 6916 K/T) cm 3 mol 1 s 1 (9) 

with uncertainties of ± 25%. The values of k<i from the present study are 
smaller than Eq. (9) by ca. 13% in the temperature interval 1400 - 1900 K, but 
they agree with each other within the combined error bounds. Because 
calculation of the temperature at [0],g for each experiment was not possible 
from the information given in Ref. 6, the source(s) of these discrepancies could 
not be found. 

Pamidimukkala and Skinner [7] determined &2 by measuring [O] with 
atomic resonance absorption spectroscopy in reflected shock heated 
H2/N20/Ar mix tures over the temperature range 1919 - 2781 K. The 
expression reported is given by 


&2 = 2.3 x 10 14 exp(- 6916 K/T) cm 3 mol 1 s l . 


( 10 ) 
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Since Eqns. (9) and (10) are basically identical, although the experimental 
temperature range is different, they are plotted as a single line in Figure II. 6. 
In the data analysis for deriving k 2 values, the converted O-atom 
concentration profiles from absorption traces with the pre-determined 
calibration curve were simulated using the reaction mechanism which was 
composed of 11 reactions. In the temperature range 1900 - 2400 K, Eq.(10) 
yielded about 20% higher values than the present results. 

Frank and Just [25] measured k 2 using atomic resonance absorption 
spectroscopy to monitor H-atoms formed in the shock heated N20/H2/Ar 
mixtures and suggested that the k 2 values of Pamidimukkala and Skinner 
should be reduced by 20%. The k 2 expression reported by Frank and Just is 

k 2 = 1 -85 x 10 14 exp(- 6976 K / T) cm 3 mol 1 s" 1 . (11) 

The present results (Eq.(8)) agree with Eq.(ll) within our experimental 
uncertainties in the common temperature range, Eq.(ll) being smaller by 8% 
at 1700 K and 6% at 2400 K. 


Sutherland et al. [8] determined k 2 by two different experimental 
methods: the flash photolysis-shock tube technique(FP-ST) combined with 
atomic resonance absorption spectroscopy over the temperature range 880 K 
to 2495 K, and the flash photolysis - resonance fluorescence (FP-RF) technique 
from 504 K to 923 K. O-atoms were produced by the flash photolysis of NO in 
the pseudo-first-order condition([0] « [H 2 D- The values of k 2 were determined 
from subsequent O-atom decay through the relationship of the Beer-Lambert 
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law, i.e., ln{ln{ I/I 0 » = - kist t + C and ki st = k 2 [H^ . The results of the FP-ST 
work were presented in two-parameter Arrhenius form, 

k 2 = (1.87 ± 0.12) x 10 14 exp(- 6854 ± 85 K/ T) cm 3 mol -1 s -1 . (12) 

For the FP-RF work, the Arrhenius expression is given by 

k 2 = (4.34 ± 0.25) x 10 13 exp(- 5249 ± 36 K/ T) cm 3 mol -1 s _1 . (13) 

The errors are given as one standard deviation. Examining Eq.(12) and Eq.(13) 
reveals that the 0( 3 P) + H 2 reaction shows Arrhenius behavior for each 
temperature range. The two expressions above were then combined with the 
low temperature results (297 K < T < 47 1 K) of Presser and Gordon [13] and 
presented with a three-parameter non-Arrhenius fit, 

k 2 = 5.06 x 10 4 T 2 67 exp(- 3165 K/ T) cm 3 mol" 1 s' 1 (14) 

with the estimated error of about ± 30% over the entire temperature range 
from 297 to 2495 K. In the temperature range of the FP-ST work, the non- 
Arrhenius expression (Eq.(14)) gives 30% higher value at 880 K and 40% at 
2495 K than Eq.(13) (Figure II.7). This increasing discrepancy toward the high 
temperatures indicates that Eq.(14) overemphasized the low temperature 
results of Presser and Gordon. As seen in Figures II. 6 and 7, the results of the 
present study agree remarkably well with Eq.(12), Sutherland et al.’s high 
temperature shock-tube results. 
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Figure II. 7 Comparison of the present rate coefficient fit of with the 
experimental data, high and extended temperature fits of Sutherland et al. 
[Ref. 8] . O is the experimental data of Sutherland et al. 
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Natarajan and Roth [9] studied reaction (2) behind reflected shock 
waves in the temperature range of 17 13 - 3532 K at total pressures of about 
1.4 - 2.0 atm by atomic resonance absorption spectroscopy using N 20 /H 2 /Ar 
mixtures. The time dependent O- and H-atom concentrations were monitored. 
A reaction mechanism of 13 reactions with appropriate rate coefficient 
expressions was used in the simulation of the experimental profiles. Through 
sensitivity and reaction flux analyses, thermal decomposition of N 2 O and 
reaction (2) were identified as the most important reactions and accordingly 
those two rate coefficients were varied to match the experimental profiles. The 
rate coefficient of N 2 O thermal decomposition, &(N20 + M), was varied within 
the uncertainty limits of ± 20%. The results were presented by a non- 

Arrhenius expression, 

k 2 = 3.72 x 10 6 T 2 17 exp(- 4080 K/T) cm 3 mol" 1 s -1 ( 15) 

with uncertainties of ± 9%. Although the values of Eq.(15) agree with the 
present results within the mutual error limits at T < 2100 K, Eq.(15), likewise 
Eq.(14), gives higher values towards the high temperatures. For example, at 
2400 K, it is 37% higher. The discrepancy between our results (Eq.(8)) and 
Sutherland et al.'s extended temperature fit (Eq.(14)) was resulted from pure 
fitting error whereas the 37% discrepancy between Eq.(8) and Eq.(15) is real. 
We repeated simulations for the maximum O-atom concentration, [0] max , and 
the time to reach [O] mox , t max , using their reaction mechanism and rate 
coefficient expressions for conditions at 1713 K (Figure 2(a)), 2171 K (Figure 
2(b)), and 2719 K of Ref. 9 where W and [O]^ are listed. The value of k 2 
was fixed to that reported by Natarajan and Roth and &(N20 + M) was varied 
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to match [O]^* exactly and compared t max . For their Figure 2(a) and (b) 
conditions, about 8% increase of &(N 2 0 + M) resulted in reasonable matches to 
the tmax values by considering ± 10 (is error on time zero of their experiments. 
However, at 2719 K, it was necessary to reduce &(N 2 0 + M) by 35%, which is 
clearly outside the error limits given for fc(N 2 0 + M). When we matched t max 
values exactly and examined the resultant maximum O-atom concentrations, 
the agreement was not much improved. Substituting their reaction 
mechanism with our Table II. 3 mechanism did not improve the matches. The 
only difference in the simulations is the thermodynamic properties of the 
species involved. Because the thermodynamic properties of the species used 
are well known, there would be no difference in computing the shock 
parameters and equilibrium constants. Frank and Just [25] used exactly the 
same method as Natarajan and Roth and reported even lower values of & 2 
than those of the present study and Sutherland et al. Therefore, the 
discrepancies at high temperatures may be due to the calibration error for the 
absorption cross sections of 0- and H-atoms used in extracting [H] and [O] in 
the experiment of Natarajan and Roth. 

Shin et al. [10] determined k<i by different experimental technique than 
atomic resonance absorption spectroscopy for O- or H-atoms. A laser 
absorption spectroscopic method was used to detect the OH radicals generated 
by reflected shock heating of lean H 2 /0 2 /Ar mixtures. The maximum slope 
measured from each experimental absorption record was converted to the 
maYimnm growth rate of OH, (d[OH]/dt W, with the absorption coefficient of 
OH, which was obtained by self-calibrating the observed signal at partial 
equilibrium and the partial equilibrium OH concentration( [O H] pe) computed 
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with a reaction mechanism of 9 reactions. In the computer simulation, k 2 
values were adjusted to match the experimental (d[OH]/dt) ma ^- The reported 
k 2 expression is 

k 2 = 7.9 x 10 14 exp(- 9382 K / T) cm 3 mol" 1 s" 1 ( 16) 

with the data scatter of ± 15% over the temperature range 1790 to 2250 K. 
Although Eq.(16) agrees with the present study within the combined error 
limits (± 27%) at temperatures below 2100 K, it gives about a 40% higher 
value at 2250 K. The partial equilibrium concentration of OH computed from 
the reaction mechanism of Shin et al. is close to that of the Table II. 3 
mechanism. The maximum OH growth rate should be independent of possible 
contaminants and of the absolute value of the absorption coefficient used in 
the simulation and experimental data evaluation because the self-calibration 
scheme was applied to both cases. We checked whether the high k 2 values of 
Shin et al. came from the reaction rate coefficients used in their reaction 
mechanism in the following way: An experimental condition, T=2250 K and 
P=1.4 atm, was chosen arbitrarily and (dfOHJ/dt),^ and A mar (instead of Ape) 
were calculated with their reaction mechanism. These two values were 
assumed to represent the experimental data. The rate coefficients k\, k%, and 
&4 were then replaced by those in the Table II. 3 mechanism and the 
(d[OH]/dt) ma * was recalculated while keeping A max the same as before. The 
maximum OH growth rate decreased by only 1.4%. This indicated that the 
difference between the reaction mechanism of Shin et al. and that of the 
present study is the k 2 expression itself. Unfortunately, we were not able to 
find the source(s) of these discrepancies because of the lack of information 



100 


given in their paper. 

Recently, Yang et al. [26] repeated the experiments to determine k 2 
values using the improved OH laser absorption technique used by Shin et al 
[10] behind the reflected shock waves; this time the laser wave-length was 
fixed by an intracavity assembly (ICA). Besides the same H 2 /O 2 mixture as 
used before (1750 K - 2250 K), a 0.05% N 2 0, 0.1% H 2 , and 99.85% Ar mixture 
(1600 K - 2140 K) was also used. Again the maximum OH growth rate was 
modeled using the reaction mechanism of 14 reactions, including N 2 0 
reactions. The new reaction rate coefficient expression for k 2 reported was 

k 2 = 3.63 x 10 14 exp(- 7818 K / T) cm 3 mol -1 s' 1 ( 17) 

with uncertainties of ± 15% from 1600 to 2250 K. Eq.(17) agrees with Eq.(8) 
within the combined error limits, Eq.(17) being larger than Eq.(8) by 6% at 
1600 K and by 26% at 2250 K. However, as shown in their sensitivity analysis 
for experiments with the N 2 0/H 2 /Ar mixture, &(N 2 0 + M) has the highest 
sensitivity on (d[OH]/dtW (Figure 2(b) in Ref. 26). Therefore, we examined 
the effect of £(N 2 0 + M) on the evaluation of k 2 in the following manner. We 
as s um ed the Table 1 reaction mechanism of Yang et al., the OH absorption 
coefficient of Shin et al.(Ref. 10) correctly regenerated the experimental 
(d[OH]/dt) max and A max at the condition of Figure 2(b) of Ref. 26 (1850 K, 1.53 
atm). Then the k 2 expression of Eq.(17) was replaced by our results (Eq.(8)) 
and &(N 2 0 + M) and e(OH) were varied to match (dtOHl/dfW* and A majc . Little 
change in e(OH) was necessary but &(N 2 0 + M) had to be increased by 10%. 
The k(N 2 0 + M) used in the Table 1 mechanism of Ref. 26 is about 7% small at 
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1850 K compared to that of Frank and Just [25], which had a ± 20% 
uncertainty limit. In other words, if we use Frank and Just's &(N20 + M) and 
increase it by only 3%, our &2 value also predicts the experimental 
observations at this temperature quite well. Furthermore, reaction flux 
analyses show that about 40% of the OH formed at the time of (d[OH]/dt) mai is 
from the reaction of N 2 O + H = N 2 + OH and the rate coefficient of this reaction 
is unknown at high temperatures (See Figure 16, p.404 of Ref. 27). At 2140 K, 
Eq.(17) is 24% higher than Eq.(8) and it was necessary to raise k(N 2 0 + M) of 
Frank and Just by 14% to obtain a satisfactory match. A 14% variation of 
fe(N20 + M) is still within its error limit, but seems too large. Obviously, high 
experimental (d[OH]/dt) max caused a high £2 value. Again, it was not possible 
to determine the cause of this difference because the experimental conditions 
were not given for their oscillogram. 

Davidson and Hanson [11] performed experiments to measure the rate 
coefficient of reaction (2). Mixtures of N0/N20/H2/Ar were first shock-heated 
and then photolyzed with an ArF excimer laser. The pyrolysis of N 2 O and the 
photolysis of NO were utilized to produce O-atoms under near identical 
conditions. The measured O-atom profiles by atomic resonance absorption 
spectroscopy were used to derive the k 2 values. Since the experiment with O- 
atom generation by NO photolysis was done under pseudo-first-order 
conditions, the slope of ln{ln{ I/I 0 ))t was equated to £ 2 ^ 2 ] and subsequently k 2 
was evaluated with known [H 2 L The Arrhenius expression for k 2 hi the 
temperature range 2120 to 2750 Kis 


k 2 = 8.13 x 10 14 exp(- 9540 ± 800 K/T) cm 3 mol 1 s 1 


(18) 
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with the uncertainties in the A factor of ± 10%. The rate coefficients 
determined by the two methods were in essence indistinguishable. Compared 
to Eq. (8), the k 2 values from Eq. (18) are about 35% high in the overlapping 
temperature range. Because Eq. (18) has higher activation energy than Eq. 
(8), the value of Eq. (18) is about 70% higher than that of Eq. (8), extrapolated 
to 2750 K. As pointed out in the literature [8], the k 2 values at high 
temperatures should be higher than those obtained by extrapolation of the low 
temperature expression because of the contribution from the reaction of 
0( 3 P)+H 2 (i; > 1). To find the possible reason(s) of discrepancy, computer 
simulations were performed under the condition of Figure 1 of Davidson and 
Hanson [11]. The reaction mechanism of Masten et al. [22] was modified by 
replacing the k 2 expression by Eq.(18) and adding N/O, N/H, and H/N/O 
reactions with rate coefficient expressions from Hanson and Salimian [27] . In 
the simulation, c(O) of 14000 atm' 1 cm' 1 and absorbance of NO of 0.158 at 
130.5 nm were used as provided [28] to convert the absorption profiles to O- 
atom concentration profiles. 

For the photolysis experiment, Eq. (18) predicts O-atom decay quite well 
while Eq. (8) generates more O-atoms with a slightly slower decay rate. During 
the course of simulations it was found that [0] ma * and t ma x (time to reach 
[0] max ) were dependent upon the reaction rate of N + NO = N 2 + O. When the 
rate coefficient of Davidson and Hanson [29] for this reaction was used, a 
much shorter t max with a higher [0] max (about 20%) than the experimental 
observations resulted. For the pyrolysis experiment, if N 2 0 was used as an O- 
atom source, [0] max was too low and t max was too long compared to the 
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experiment. In order to match the experimental O-atom profile, first, some 
other O-atom source, e.g., NO 2 instead of N 2 O should be used as an O-atom 
source [28]; second, the initial concentration of NO 2 should be higher than the 
N 2 O concentration appearing in the caption of Figure 1 in Ref. 11. Accordingly, 
the initial NO 2 concentration was adjusted to give an excellent match to the 
converted experimental O-atom profile. After this was done Eq.(18) was 
replaced by Eq.(8). Again, with Eq.(8), high O-atom concentrations with a 
slightly slower decay rate compared to the experiment could be observed. 
Since the mixture compositions are not known exactly, it was difficult to 
conclude that the k<i expressions obtained from pyrolysis and photolysis 
methods were indistinguishable. For h<i evaluation, the most important 
element is the absorption cross section of O-atom. It has been well known that 
calibration for the atomic species using the atomic resonance atomic 
spectroscopy method is extremely difficult. For example, the absorption cross 
sections of H(D)- and O-atom vary depending upon the construction and run 
condition of microwave discharge lamps [30, 31, 32, 33]. Therefore, it may be 
better to use the self-calibration scheme to match the experimental absorption 
profiles. If g(O) is reduced by about 35% from the value used by Davidson and 
Hanson, then our &2 expression (Eq.(8)) would also give an excellent match to 
the experimental profiles. Furthermore in Davidson and Hanson's 
experiments, especially at high temperatures, the flow effect could not be 
ignored because NO was photolyzed at long time after the shock wave passage 
(ca. 550 ps) and subsequent decay of O-atoms was used for the evaluation of 
&2 values. In such cases, a possible change of flow conditions should be 
properly considered. There is also a possibility of O-atom depletion to the 
boundary layer which causes a fast O-atom decay because very small O-atom 
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concentrations were monitored throughout the experiments. Nonetheless, 
even though considering the enhancement of the rate coefficient at high 
temperatures by vibrational excitation of H 2 , our rate coefficient expression 
yields considerably smaller values as compared to Eq.(18) at high 
temperatures. 

One of the most thorough theoretical estimates made on reaction rates 
is for reaction of 0( 3 P) + H 2 = OH + H. Bowman et al. [34] calculated this 
reaction rate coefficient using the reduced-dimensionality collinear exact 
quantum with bend (CEQB/G) reaction probability and transmission 
coefficient based on an effective potential energy surface given by the ab initio 
MOD POLCI potential energy surface [35] plus ground state adiabatic bending 
energy. The calculated rate coefficients were fitted to an Arrhenius expression 
in the temperature of 1400 to 1900 K [12] and were given by 

k 2 = 2.89 x 10 14 exp(- 7240 K/ T) cm 3 mol" 1 s" 1 ( 19) 

Equation (19) gives about 16% and 26% higher values than Eq.(8) at 
temperature of 1400 K and 1900 K, respectively (Figure II.8). Previous 
calculations by Bowman and his coworkers using the collinear exact quantum 
(CEQ) reaction probability and transmission coefficient based on MOD POLCI 
surface without incorporating ground state bending energy [35] yielded slightly 
lower values than the CEQB/G method [34], and hence gave better matches to 
the experiments. Recently Truhlar and coworkers also calculated the rate 
coefficients for the 0( 3 P) + H 2 = OH + H reaction using the ICVT/LAG 
(Improved Canonical Variational Transition State Theory with Least-Action 
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10kK/T 


Figure II. 8 Comparison of the present results for k<i with the high 
temperature and the extended temperature fits of Sutherland et al.[Refi 8] 
and the theoretical calculation of Bowman et al.fRef. 34, 35] and of Truhlar 
and his coworkers [Ref. 12, 36] . 




106 


Ground-state transmission coefficient) method [36 and references therein]. 
The potential energy surface used (J3 surface) was the modification of RMOS- 
MOD POLCI surface for collinear O-H-H geometries of Lee et al. [35] 
augmented by including bend potentials for 3 A" and 3 A' surfaces. The classical 
barrier height was recalibrated to 13.0 kcal/mol to match to the experiments. 
The rate coefficients computed at high temperatures (in units of cm 3 mol ^s x ) 
are 2.53 x 10 11 , 2.11 x 1()12, and 1.33 x 10 13 at 1000 K, 1500 K and 2400 K, 
respectively (Figure II. 8). At low temperatures (T < 500 K), there were 
excellent matches to the experiment of Presser and Gordon [13] within 7%, the 
computation being 7% higher at 472 K. Inspection of the results indicates that 
if a curve is used to fit the data, it is highly curved upward. Compared to the 
results of the present study, computations give about 8% higher values at 
1500 K and 24% higher at 2400 K. The discrepancy between the present 
results and the calculations is possibly due to the assumption of non-recrossing 
of trajectories and/or too much reaction flux through the 3 A' potential energy 
surface at high temperatures. 


Summary 

The measurement of the rate coefficients for 0( 3 P) + H2 = OH + H 
reaction by OH laser absorption spectroscopy coupled to a shock-tube 
technique has been described and discussed. 


In the experiments, actual pressures behind the reflected shock waves 
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were measured and used together with the adiabatic equation of state to 
correct the calculated ideal shock properties. The corrected temperatures 
were about 1.4% higher than the ideal shock temperatures. Hence this 
temperature correction has the effect of lowering the reaction rate coefficients. 
Experimental conditions were chosen such that the 0( 3 P) + H 2 = OH + H 
reaction gave the highest sensitivity to the experimental observables used in 
the optimization of the rate coefficient. The OH absorption profiles were 
successfully reproduced with Table II. 3 reaction mechanism containing k -2 
expression below. 

The present data can be best fitted by an Arrhenius expression in the 
temperature range 1424 K to 2427 K: 

k 2 = (1.88 ± 0.07) x 10 14 exp(- 6897 ± 53 K/ T) cm 3 mol" 1 s" 1 (20) 

with the error limits of ± 12%. The results of this study agree quite well with 
the high temperature rate coefficient expression of Sutherland et al. [8]. 
However, Eq.(8) disagrees with the measurement of Natarajan and Roth, Shin 
et al., Davidson and Hanson, and Yang et al. who used N 2 O as an O-atom 
source [9, 10, 11, 26] and simulated either the O-atom concentration profile 
converted from the absorption profile with the absorption cross section of O- 
atom [9, 11] or maximum OH formation rate obtained from the absorption 
profile with the OH absorption coefficient [10, 26]. The discrepancies can be 
resolved or lessened if one considers the difficulty in calibrating the O-atom 
absorption cross section and the uncertainty of the N 2 O decomposition rate. It 
is worthy to note that the extended temperature expression of Sutherland et al. 
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[8] gives considerable overprediction of k 2 values at high temperatures 
compared to both our and their own high temperature data. This suggests that 
the upward curvature of the &2 expression due to tunneling at low 
temperatures and contribution from 0( 3 P) + H2 (u> 1) at high temperatures 
may not be adequately represented by a 3-parameter non- Arrhenius equation. 

The theoretical calculations [12, 34, 361 are in good agreement with the 
extended temperature expression of Sutherland et al., which is the combined 
results of the low temperature data of Presser and Gordon [13] and middle to 
high temperature data of Sutherland et al. [8]. Since the extended 
temperature expression of Sutherland et al. overpredicts the high temperature 
data, theoretical calculations also overestimate the high temperature reaction 
rates. Further theoretical developments would be required to resolve the 
existing discrepancies between the experimental and calculated ^2 values 
especially at high temperatures. 
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